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ABSTRACT 

' We report on a recent event in which, after more than a decade of slowly fading, the 

visual brightness of the massive young binary Z CMa suddenly started to rise by about 
CO . 1 magnitude in December 1999, followed by a rapid decline to its previous brightness 

over the next six months. This behaviour is similar to that exhibited by this system 
around its eruption in February 1987. A comparison of the intrinsic luminosities of 
the system with recent evolutionary calculations shows that Z CMa may consist of 
a 16 M Q BO Hie primary star and a ~3 M© FUor secondary with a common age 
of ~ 3 x 10 5 yr. We also compare new high-resolution spectra obtained in Jan. and 

■ Feb. 2000, during the recent rise in brightness, with archive data from 1991 and 1996. 

The spectra are rich in emission lines, which originate from the envelope of the early 

q . B-type primary star. The strength of these emission lines increased strongly with 

!— i ' the brightness of Z CMa. We interpret the collected spectral data in terms of an 

\ accretion disc with atmosphere around the Herbig BOe component of Z CMa, which 

■ has expanded during the outbursts of 1987 and 2000. A high resolution profile of the 

6300 A [O i] emission line, obtained by us in March 2002 shows an increase in flux and 
a prominent blue shoulder to the feature extending to ~ —700 km s , which was 

■ much fainter in the pre-outburst spectra. We propose that this change in profile is a 

result of a strong change in the collimation of a jet, as a result of the outburst at the 
start of this century. 

Key words: circumstellar matter - stars: emission-line - stars: pre-main sequence - 
stars: Z CMa - stars: variables - open clusters and associations: CMa Rl 



o 
o 



* Based on observations made with the European Southern Ob- 
servatory telescopes obtained from the ESO/ST-ECF Science 
Archive Facility, with the William Herschel Telescope operated 
on the island of La Palma by the Isaac Newton Group in the 
Spanish Observatorio del Roque de los Muchachos of the Insti- 
tuto de Astrofi'sica de Canarias, with the BTA telescope of the 
Special Astrophysical Observatory, Nizhnii-Arkhyz, Russia, and 
with the Alfred Jensch telescope of the Thiiringer Landesstern- 
warte, Germany. 

' Visiting Astronomer at the Infrared Telescope Facility, which is 



1 INTRODUCTION 

Ever since the discovery of bright emission lines in its 
spectrum by Merrill (1927), Z CMa (HD 53179, IRAS 
07013—1128) has been one of the most frequently studied 
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ment no. NCC 5-538 with the National Aeronautics and Space 
Administration, Office of Space Science, Planetary Astronomy 
Program. 
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massive young stars. The brightness of the system is highly 
variable, varying irregularly between its quiescent state at 
V « 11 to its active state at V w 9 (Covino et al. 1984; 
Hessman et al. 1991). Although the photometric and spec- 
troscopic properties of Z CMa show some similarity with 
those of the FUors (e.g. FU Ori and V1057 Cyg), its spec- 
tral differences with these objects, especially in the presence 
of many narrow emission lines, are remarkable. Nevertheless, 
in its normal brightness phases (V > 9.6 m ) the absorption 
character of the spectrum is dominant and, similar to the 
FUors, the spectrum can be interpreted as that of an opti- 
cally thick accretion disc (Hartmann et al. 1989; Welty et 
al. 1992). 

In Feb. 1987 a large rise in brightness of Z CMa was ac- 
companied by a rise in emission flux in the Balmer lines and 
CallK and with the transition of many metal lines of Fen, 
Crn and Till from absorption into emission (Hessman et 
al. 1991). These authors tried to explain the absorption and 
emission spectral components in terms of a common source, 
a variable accretion disc. However, it proved difficult to find 
a simple model for such a source that would be able to re- 
produce the observations. A key ingredient in solving this 
mystery was discovered in 1989, when near-infrared speckle 
interferometry revealed that Z CMa is a double star with 
a separation of 0.l"(Koresko et al. 1991; Haas et al. 1993). 
These authors showed that the South-Eastern component of 
Z CMa dominates the visual and UV parts of the SED and 
is associated with a FUor-like system, whereas the North- 
western component is a powerful infrared source, which is 
the primary object of this double star system. The separa- 
tion of Z CMa into two components has since been confirmed 
by speckle masking observations in the visual (_R-band) by 
Barth et al. (1994) and by Thiebaut et al. (1995). In addition 
it was discovered by Whitney et al. (1993) that the emission 
lines in the visual part of the spectrum of Nov. 1991 were po- 
larized, in contrast to the continuum of the spectrum. This 
supports the idea that these emission lines are contributed 
by scattered light from the primary source of the system. 

The analogy of this spectrum with that of the early 
type Herbig Be star MWC 1080 suggests that the source 
of the emission lines is an early type Be star inside a large 
cocoon. The radiation of this Be star escapes by scattering 
from the dust walls of a cavity in the thick cocoon envelope. 
Millan-Gabet & Monnier (2002) succesfully imaged the di- 
rect circumstellar environment of Z CMa at 1.25 fxm, and 
found a small (~ 1") linear feature. They interpreted this 
feature as light scattered off the walls of a jet-blown cavity, 
confirming the above picture. The cavity could be formed by 
the extended (over 3.6 pc) optical jet discovered by Poetzel 
et al. (1989). This jet is also seen as a bipolar outflow in CO 
(Evans et al. 1994). More recent spatially resolved optical 
spectroscopy by Bailey (1998) and Garcia et al. (1999) has 
confirmed that the extended outflow originates from the in- 
frared primary component of Z CMa and that the optical 
emission line spectrum is associated with the same compo- 
nent. 

During Nov. /Dec. 1999, observers of the ROTOR pho- 
tometric monitoring program at Mt. Maidanak Observatory 
(Uzbekistan) found that after seven years of low brightness 
(V ~ 10. 3 m ) the brightness of Z CMa rose sharply (within 
40 days) by more than 1 magnitude. We obtained new high- 
resolution spectra of Z CMa shortly after this new outburst. 



These spectra show intrinsic changes in line profiles and 
equivalent widths with respect to spectra obtained in 1991 
and 1996. In this paper we discuss the spectral changes of 
Z CMa over the last ten years and interpret these changes 
as due to changes in the size of an accretion disk, associated 
with the Herbig BOe primary. 

From UBVRI photometry of Z CMa obtained in Nov. 
1991 and Feb. 1987 (Sect. 3.2) we confirm the spectral clas- 
sification of the star in the primary source and derive its 
radius and luminosity for spectral type B0III. Although the 
photometric data of Z CMa suggest that the variability of 
the contributions from both binary components is for a large 
part due to irregular variations in their circumstellar dust 
extinctions (Sect. 3.1), there are variations in the spectrum 
of the Be component which show that intrinsic changes of 
the Be envelope excitation occur and possibly also influence 
the circumstellar extinction. 

In Sect. 4 we discuss the formation regions of the vari- 
ous emission lines, and note that several strong absorption 
lines from He I, Ol and part of the NalD lines can also be 
attributed to the BOe star. From the red Ca II (2) emission 
line triplet and various multiplets of Fe II emission lines we 
have estimated lower limits to the disc radius of the 'hidden' 
Herbig Be star. We also make estimates of the outflow and 
accretion rates of the Be star. 

In the discussion (Sect. 5) we compare our data of 
Z CMa with those of three Herbig stars with similar spectra: 
MWC 1080, V645 Cyg and V380 Ori. We also make some 
attempts to estimate the evolutionary stage and mass of the 
BOe star from recent models for the pre-main sequence evo- 
lution of massive stars. It appears possible to find models for 
the two components of Z CMa which predict the accretion 
rates (derived from the observations) and which also predict 
a bipolar jet only for the BOe star. 



2 OBSERVATIONS 
2.1 Photometry 

Since 1980 Z CMa has been monitored closely in the John- 
son UBVR system by the ROTOR photometric monitor- 
ing program. The observations in this program were made 
with three different 0.6 m telescopes at Mt. Maidanak Ob- 
servatory, Uzbekistan. During these observations, mostly 
made through a 15" circular diaphragm, the telescope was 
equipped with a single-channel photometer. Standards and 
extinction reference stars were observed each night. The 
data were reduced with standard techniques at the Tashkent 
Astronomical Institute. A detailed description of the RO- 
TOR program and the employed measuring procedures can 
be found in Shevchenko (1989). 

The ROTOR data on Z CMa were supplemented with 
optical and infrared measurements on various photometric 
systems from the literature. When possible, all data were 
transformed to the "standard" Johnson/ Cousins UBVRcIc 
and near-infrared JHKLM systems using the transforma- 
tion formulae from Fernie (1983), Brand & Wouterloot 
(1988), and Carpenter (2001). The results are given and dis- 
cussed in Sect. 3.1. 
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Table 1. Log of spectroscopic observations of Z CMa. We also list the V-band photometric measurement closest to the date of each 
spectroscopic observation. 



Telescope Range [A] Date JD+2400000 R [mA] t oxp [min] V [mag] 



SAO 6 m 


4190-5490 


22 


Nov. 


1991 


48582.83 


300 




9.7 


SAO 6 m 


5960-7150 


23 


Nov. 


1991 


48583.83 


300 




9.7 


ESO NTT 


5720-9930 


20 


Dec. 


1991 


48610.838 


279.5 


20 


9.8 


ESO NTT 


5720-10030 


21 


Dec. 


1991 


48611.710 


279.5 


45 


9.8 


ESO 1.5 m 


4200-7800 


6 


Dec. 


1992 


48962.823 


1076 


4 


10.3 


ESO 1.5 m 


6200-8700 


10 


Dec. 


1992 


48966.651 


747 


8 


10.3 


ESO CAT 


5860-5910 


15 


Dec. 


1994 


49701.355 


107.0 


45 


10.2 


ESO CAT 


6521-6600 


12 


Dec. 


1996 


50429.151 


65.6 


45 


10.2 


ESO CAT 


5848-5918 


13 


Dec. 


1996 


50430.210 


58.8 


45 


10.2 


ESO CAT 


5848-5918 


14 


Dec. 


1996 


50431.235 


58.8 


45 


10.2 


ESO CAT 


3913-3955 


14 


Dec. 


1996 


50431.323 


39.3 


45 


10.2 


ESO CAT 


8467-8573 


15 


Dec. 


1996 


50432.235 


85.2 


45 


10.2 


WHT 


5220-9110 


26 


Dec. 


1996 


50443.725 


266.7 


20 


10.3 


WHT 


5220-9110 


26 


Dec. 


1996 


50443.740 


266.7 


20 


10.3 


Tautenburg 2 m 


5600-10050 


23 


Jan. 


2000 


51567.422 


230.0 


20 


9.1 


WHT 


3640-5470 


18 


Feb. 


2000 


51593.361 


177.2 


30 


9.2 


WHT 


4100-8350 


18 


Feb. 


2000 


51593.423 


232.0 


30 


9.2 


WHT 


4600-9960 


18 


Feb. 


2000 


51593.447 


278.4 


30 


9.2 


WHT 


5330-11020 


18 


Feb. 


2000 


51593.498 


325.8 


30 


9.2 


NASA IRTF 


19000-41400 


18 


Dec. 


2001 


52262.139 


32.0 


5 


10.0 


ESO 3.6 m 


6279-6321 


31 


Mar. 


2002 


52365.084 


42.0 


30 


10.0 



2.2 Spectroscopy 

In this paper we present the spectra of Z CMa obtained with 
a range of telescopes: 

(a) Two low resolution spectra, taken on 9 and 10 De- 
cember 1992 with the Boiler and Chivens spectrograph on 
the 1.5 m telescope at ESO's La Silla Observatory. These 
spectra show emission in Ha, [Ol] 6300 A, O I 8446 A and 
Can (2) (8500, 8545 and 8662 A). The emissions of H/3, P14 
and Fell (42), which are clear in the low resolution spectrum 
of the B5 star LkHa 220 (Tjin A Djie et al. 2001 (here- 
after paper II), Fig. 5) are much weaker in this spectrum of 
Z CMa. In contrast, the emission lines in the Z CMa spec- 
tra of Whitney et al. (1993, Fig. 1). of Nov. 27, 1991 are 
stronger than those of LkHa 220. 

(b) A high resolution blue spectrum (4190-5490 A) ob- 
tained with the Main Spectrograph on the 6 m BTA tele- 
scope of the SAO (Special Astrophysical Observatory in 
Nizhnii-Arkhyz) on 22 Nov. 1991 and a high resolution Ha 
profile obtained with the same instrument on 23 Nov. 1991. 
Similar to the spectrum of Whitney et al. (1993), which was 
taken a few days later, the SAO spectrum shows strong emis- 
sion components in Ha, H/3, H7 and in several multiplets of 
Fen. These spectra are supplemented by high resolution red 
ESO Multiple-Mode Instrument (EMMI) spectra, obtained 
on 20 and 21 Dec. 1991 with the 3.5 m ESO New Technology 
Telescope (NTT) and extracted by us from the ESO science 
archive. 

(c) High resolution echelle spectra have been secured 
with the Coude echelle spectrograph (CES) mounted on the 
1.4 m Coude Auxiliary Telescope (CAT) at La Silla: (1) 
NalD absorption profiles of Dec. 15, 1994 and of Dec. 13 
and 14, 1996 (2) An Ha emission profile of Dec. 12, 1996 (3) 
A CallK emission and absorption profile of Dec. 14 1996 
(4) Can (8498 A) and Call (8542 A) (blended with resp. 
P16 and P15) emission profiles of Dec. 15, 1996. Addition- 



ally, an emission profile of the [0 1] 6300 A line has been 
obtained with the same instrument connected to the ESO 
3.6 m telescope on March 31, 2002. 

(d) A high resolution red spectrum (5470-9110 A) taken 
on Dec. 26, 1996 (close to the dates of the observations men- 
tioned under c) with the Utrecht Echelle Spectrograph (UES) 
on the 4.2 m WHT at the Northern Hemisphere Observa- 
tory on La Palma was kindly made available to us by Drs. 
H. van Winckel and G. Meeus of the Katholieke Universiteit 
Leuven (Belgium). 

(e) A high resolution red spectrum (5600-10050 A) 
was obtained on January 23, 2000 with the Coude Echelle 
Spectrograph of the 2 m Alfred Jensch telescope of the 
Thiiringer Landessternwarte 'Karl Schwarzschild' in Taut- 
enburg, Thiiringen. This spectrum was taken during a very 
bright phase of Z CMa (V ~ 9.1 m ), comparable to the out- 
burst of Feb. 1987 (V ~ 8.7 m ). The spectrum covers the 
wavelength range from about 560 to 1005 nm. Using a Tek 
CCD Chip with 1024 x 1024 pixels and a slit-width of two 
arcseconds, the two-pixel resolution of the spectrum is about 
A/ (<5A) = 35,000. Standard IRAF routines were used in order 
to extract and wavelength calibrate the spectrum. 

(f) Four high resolution echelle spectra, centered at 
4300 A, 5200 A, 6200 A and 8200 A were obtained on Feb. 
18, 2000 with the UES on the 4.2 m WHT. A number of line 
profiles from these spectra will be compared with the corre- 
sponding profiles from the spectra mentioned under (b)-(e) 
and with those published in the literature. Apart from the 
blue range, these spectra extend in the red up to 11020 A. 
However, some lines (e.g. [0 1] 6300 A) could not be observed 
in these spectra because there is no wavelength overlap of 
the successive orders. 

(g) A near-infrared (1.9-4.1 /jm) spectrum has been 
taken with the SpeX spectrograph (Rayner et al. 2003) 
attached to the 3.0 m NASA Infrared Telescope Facility 
(IRTF) of the Mauna Kea Observatory on Dec. 18, 2001. 
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Figure 1. Optical (UBVRcIc) light-curve of Z CMa spanning the 
period 1980-2003. Typical errors in magnitude and colours are of 
the order of 0.02 mag. The triangles in the top panel indicate the 
dates of the spectroscopic observations listed in Table 1. 



In this spectrum we observe Br a, Br 7, Pf 7 and several 
CO lines in emission. 

Observational details of all spectroscopic observations 
are reviewed in Table 1. 



3 PHOTOMETRY AND INTERPRETATION 
3.1 Variability and Lightcurve 

Figures 1 and 2 show the optical (UBV(RI)c) and near- 
infrared ( JHK) lightcurve of Z CMa between 1980 and 2003. 
For a detailed description of the photometric behaviour of 
Z CMa before the 1987 outburst, we refer the reader to the 
papers of Covino et al. (1984), Hessman et al. (1991), and 
Lamzin et al. (1998). From Fig. 1 we conclude that whereas 
before 1987 the system brightness appears relatively con- 
stant, from Feb. 1987 to Dec. 2000 the brightness of the star 
is steadily declining by about 0.1 m per year. Clear variability 
is seen superimposed on this trend, such as a local minimum 
around December 1992. There appears to be a strong corre- 
lation between the observed brightness and colour variations 
in Z CMa. We illustrate this point in Fig. 3, where we give 
the V vs. (U - B), (B -V), (V - R c ) and (V - I c ) dis- 
tributions for the time-interval 1980-2003. Although some 
correlation is present in these diagrams, the spread of the 
data-points is much larger than expected from the observa- 
tional errors (~ 0.02 m ). 

A closer inspection shows that between Feb. 1987 and 
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Figure 2. Near-Infrared ( JHK) light-curve of Z CMa spanning 
the period 1980-2003. Typical errors in these data are about 0.05 
mag. 



Nov. 1991 the photometric variations closely follow a be- 
haviour of redder colours when the system gets fainter, but 
that the data of Dec. 1999 and March 2001 are shifted to 
somewhat lower values of (B — V) and (U — B) than defined 
by the reddening line. We infer that either variable circum- 
stellar extinction or a slow cooling-down of the FUor com- 
ponent in Z CMa may explain the long-term overall trend 
of the photometric variation between Feb. 1987 and Dec. 
2000. However, on a shorter time-scale of weeks and days 
the variations seem erratic and hard to explain. This may 
not be surprising, since we probably are observing the re- 
sult of two independently varying sources, a FUor secondary 
component which usually dominates the visual continuum 
and a Be type primary star in a large dust cocoon, which 
is dominating the near and far infrared part of the spec- 
trum (Koresko et al. 1991). Fig. 2, displaying the variation 
of the near-infrared magnitudes, shows that in 1991 the H 
and K band fluxes reach a local maximum, which is more 
pronounced than that in the V band. This suggests that the 
circumstellar extinction of the cocoon around the primary 
component has decreased so far, that the visual flux from 
that component could increase its contribution to the total 
flux in the visual somewhat. 

Apart from the long-term variations discussed above, 
the dominant features in the optical light-curve of Z CMa are 
the large and fast rises in visual flux in Feb. 1987 and Jan. 
2000 (Fig. 1). The latter maximum started with a fast rise, 
first observed in Oct. 1999, over 0.7 m in V taking 46 days. 
The star remained at this high brightness for at least 180 
days. When extrapolating to earlier dates, the total increase 
in V-brightness is ~l.l m , but the time-scale of this larger 
variation has not been observed. The time-scale of the rise 
in brightness over 0.7 m during Jan. /Feb. 1987 is also not 
known, but the duration of the higher brightness phase was 
also about 6 months (Hessman et al. 1991, Fig. 1). 

We note that the large and fast rises in the visual fluxes 
in Feb. 1987 and Jan. 2000, do not originate in a large, 
fast decrease in circumstellar extinction; it will be shown by 
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Figure 3. Optical (UBVR c Ic) colour-magnitude diagram of Z CMa using photometric data from the period 1980-2003. Data from Jan. 
1980-Jan. 1987 arc indicated by the plusses, data from the Feb. 1987 outburst by open squares, data from Mar. 1987-Oct. 1999 by the 
crosses, data from the Nov./Dec. 1999 outburst by the open circles, and data points obtained after Feb. 2000 are indicated by the asterisk 
symbols. The arrows indicate the direction of interstellar reddening 



the spectroscopy (Sect. 4) that they are related to intrinsic 
changes of the primary Be star component. 

3.2 The SEDs in the visual 

From the spectral polarization data measured in Nov. 1991 
by Whitney et al. (1993), the authors derived that at this 
date 75% of the visual flux was contributed by the FUor 
secondary source and that the remaining 25% was due to 
the Be-type primary star. By assuming that the fast rise 
in visual brightness of Feb. 1987 (Hessman et al. 1991) was 
entirely due to a rise in the primary brightness only, they 
derived that on Feb. 19, 1987 the contribution of the primary 
Be star to the total visual flux was 66%. In the same way we 
derive from the fast rise in Nov./Dec. 1999 that at maximum 
brightness in Jan. 2000 (V ~ 9.1 m ) the contribution of the 
primary Be star to the total visual flux was ~ 60%. 

In the UV the spectrum is dominated by the FUor com- 
panion. Earlier obtained UV spectra could be matched with 
a model of an accretion disc with a late F-type spectrum 
(Hartmann et al. 1989; Kenyon et al. 1989). We have made 
a similar fit of two IUE spectra (taken on 17 April 1979 when 
V ~ 9.3 m and 19 March 1988 when V ~ 9.6 m ). Good fits 
were obtained with a F5Ib star, an 0.02 R* thick boundary 
layer with T e g — 10,500 K and an optically thick disc with 
5400 K maximum temperature. We assumed a foreground 
excess of 0.16 m with Ri s — 3.1 and a foreground CMa Rl 
cluster excess of 0.20 m with R CB = 4.2 (Shevchenko et al. 
1999; Paper I). The accretion rate for the FUor is then ~ 3 
x 10~ 5 M yr" 1 . 

From the observed UBVRI photometry, the above 
mentioned foreground excess and the relative flux contri- 
butions of a B0III primary and a FUor (F5Ib) secondary 
component we have calculated the extinction-free UBVRI 
fluxes of both components in 1987, 1991 and 2000 for various 
values of R CB (3.1, 4.2 and 6.0). We assumed here intrinsic 
colours (B - V) = -0.29 m for a B0III star and (B - V) = 



Table 2. Radii and luminosities of the BOIIIc star on 20 Feb. 
1987, 23 Jan. 2000, and 22 Nov. 1991 for R cs = 3.1, R cs = 4.2 
and R CB = 6.0. The uncertainties in the computed values of R/Rq 
and log(L/Lo), dominated by the 14% error in the distance to 
CMa Rl (paper I), are less than 10% and 3around a few thou- 
sand K. 



Obs. Date 


v obs H 


Rcs 


v H 


R/R Q 


log(L/L ) 


T«h [K] 


02/20/87 


8.8 


3.1 


5.39 


12.3 


4.81 


26,000 


01/23/00 


9.2 


3.1 


5.70 


10.6 


4.68 


26,000 


11/22/91 


9.8 


3.1 


6.14 


8.0 


4.57 


28,000 


02/20/87 


8.8 


4.2 


4.50 


18.5 


5.16 


26,000 


01/23/00 


9.2 


4.2 


4.72 


16.7 


5.07 


26,000 


11/22/91 


9.8 


4.2 


4.80 


16.1 


5.04 


28,000 


02/20/87 


8.8 


6.0 


3.03 


33.6 


5.81 


28,000 


01/23/00 


9.2 


6.0 


3.12 


32.4 


5.78 


28,000 


11/22/91 


9.8 


6.0 


2.61 


37.5 


6.01 


30,000 



0.37 m for the FUor (Schmidt-Kaler 1982). The extinction- 
free UBVRI distributions of the B-type primary star can 
be matched with the Kurucz model flux ratios of a late O or 
early B-type star, which gives us an effective temperature, 
a radius and a luminosity for each of the dates of spectro- 
scopic observations (Table 2). We will use these results in 
Sect. 4 and in the discussion on the evolutionary status of 
Z CMa (Sect. 5). 

In reality the values of R/Rq and log(L/L ) of Table 
2 may be somewhat higher because so far we have not ac- 
counted for the fact that we only receive the light from the 
star after it has been scattered and polarized by the dust 
walls of the cavity in the cocoon. The correction for this re- 
duction is uncertain, but because of the upper limit of 10% 
for the degree of polarization (Whitney et al. 1993) we es- 
timate that the scattering angle will not be larger than 20 
degrees, so that the flux reduction due to scattering will not 
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be more than 5% if the scattering grains are more or less 
aligned (see Figs. 8 and 9 of the paper of Voshchinnikov & 
Farafonov 1993). Our estimate that the scattering angle is 
small is supported by the relatively large radial velocities 
of the [N n] and [S n] emission lines in the bipolar outflow 
along the axis of the cavity (Poetzel et al. 1989). 

3.3 The mass loss of the Be primary component 

From Jan. 1981 Z CMa photometry on the Walraven 
WULBV system (de Winter et al. 2002) we can make an 
estimate of the mass loss from the early Be-type companion. 
After transformation of the flux in the Walraven V-band to 
the Johnson system we find that the visual magnitude in 
this period was 9.0 m , which suggests that the continuum of 
the B star dominated the visual part of the spectrum at that 
date. After correction for foreground extinction, the Balmer 
jump Db, derived from the Walraven photometry is 0.16 m . 
Predicted values of the Balmer jump based on model cal- 
culations by Mihalas (1966) can be found in Fig. 2 of the 
paper of Garrison (1978). For BO and log<? = 3 the theoret- 
ical value of the Balmer jump is close to 0.15 m , but for log g 
= 4 the predicted value is already 0.29 m , which implies that 
ADb for BOIII has an upper limit of 0.13 m . We estimate 
the mass loss with the method of Garrison and assume a 
temperature of 10,000 K for the outflow region. Then the 
upper limit of ADb leads to an upper limit of 1.6 x 10 37 
cm~ 6 for the volume emission measure and with an outflow 
velocity of 1000 km s _1 in the Ha emission region (see Sect. 
4) and a radius of 9.4 R Q we have an upper limit of 1 x 
10~ 4 Mq yr _1 for the mass loss in the wind. The mass loss 
is proportional with the outflow velocity in the emission re- 
gion, so that the actual upper limit can easily be twice as 
low. 



4 ANALYSIS OF THE SPECTRA 

4.1 Results 

Table 4 gives the data of the lines in the red spectra (long- 
ward 5650 A) in the Z CMa spectra for V ~ 9.1 m , for V 
~ 9.7 m and for V ~ 10.3 m . The tables contain equivalent 
widths (EWs), half- widths (FWHM) and P Cygni outflow 
velocities of the lines. In Tables 3 the corresponding data of 
the lines in the blue spectra are listed. 

In order to learn more about the structure of the en- 
velope of the early Be type primary we discuss the various 
lines and their profiles in the high resolution red and blue 
spectra of Jan. and Feb. 2002 (V ~ 9.1 m ) and then compare 
the results with the corresponding red and blue spectra, ob- 
served in Nov. (blue) and Dec. (red) 1991 (V ~ 9.6 m ), and 
with the profiles of several characteristic lines, secured in 
Jan. 1996 (V ~ 10.2 m ). 

4.2 The outburst spectra of 2000 and 1987 

4-2.1 The emission lines 

In discussing the lines we first should know of which stel- 
lar component they originate. The presence of N I (8629 
and 8683 A) emission is an indication for a hot atmosphere, 
since the lines have been only observed for (early) B stars 
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Figure 4. Hi Balmer line profiles of Z CMa. (a) Ha line profile 
obtained on 12 Dec. 1996, (b) H/3 profile obtained on 18 Feb. 2000, 
(c) H7 profile obtained on 18 Feb. 2000, (d) H<5 profile obtained 
on 18 Feb. 2000. The dashed horizontal and vertical lines indicate 
the continuum level, and the rest- wavelength of the observed line, 
respectively. 
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Figure 5. Nal and He I line profile variations in Z CMa between 
1991 and 2000. The left panel shows the 5889.95 and 5895.92 A 
NalD lines, together with the 5875.97 A He I line, whereas the 
right panel shows the spectral region around the 6678.15 A He I 
line. Strong He I absorption is only present in the Jan. and Feb. 
2000 spectra. 
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Figure 6. Ol line profile variations in Z CMa between 1991 and 
2000. The left panel shows the observed variability in the O I 
6156.41 A line, whereas the right-hand panel illustrates the be- 
haviour of the Oi 7774.17 A line. 
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(Hamann & Persson 1992b). Similarly strong, high excita- 
tion Fen 9997.57 A emission has only be observed in spec- 
tra of early Be type stars, such as 7 Cas (BOIV) (Viotti et 
al. 1998). This line is probably formed by fluorescence due 
to pumping by Lya (see e.g. Rodriguez- Ardila et al. 2002). 
Whitney et al. (1993) and Garcia et al. (1999) have shown 
that the emission lines are associated with the primary Be- 
star. The main part of this emission will come from the ex- 
tended envelope of this star, but the forbidden lines of [N 11] , 
[S 11] and [0 1] have been shown to be formed in the bipo- 
lar outflow (or a jet) from the primary (Poetzel et al. 1989; 
Garcia et al. 1999). From spectral imaging Bailey (1998) re- 
vealed that also a substantial part of the central region of 
the Ha emission is emitted by the jet. The peak of the line 
is redshifted by +82 km s _1 with respect to the laboratory 
wavelength. The wings of the line have their origin in the 
extended atmosphere of the primary Be-star. Especially the 
red wing extends very far, up to nearly 1000 km s _1 . The 
emission peak of H/3 shows a similar redshift of +75 km s _1 . 

We assume that the other Balmer (and Paschen) lines 
will be formed in the same region as Ha. A strong argument 
in favour of this is that all Balmer lines (from Ha to H15) in 
the WHT spectrum have a similar outflow velocity pattern. 
From the last resolved Balmer lines (n ~ 35) we derive from 
the Inglis- Teller relation that n e <, 2.4 x 10 11 cm~ 3 in the 
region of their formation. For the Paschen lines, which have 
strong emission components (lower continuum in red part 
of the spectrum), n — 23 which corresponds to n e w 10 13 
cm" 3 . 

The emission line of O I at 8446.5 A is expected to be 
correlated with the Ha-emission, because the 3d 3 D° state of 
Ol, which can be pumped by Ly/3 (Polidan & Peters 1976), 
cascades down succesively by the emission of 11291 A to 
3p 3 P, that of 8446 A to Ss 3 ^ and finally by that of 1302 A 
to the 3P ground state. Its EW varies from negligible at 
Apr. 1989 (V ~ 9.7 m ) to 1.6 A in Dec. 1987 (V - 9.5 m , 
Hamann & Persson 1992b), 1.2 in Dec. 1991 (V ~ 9.6 m ), 
1.5 A in Jan. 1996 (V ~ 10. 2 m , Teodorani et al. 1997) to 
1.4 A on Jan 2000 (V ~ 9.1 m ). For the last three dates, the 
EW of Ha varied from 19.1, 32.5 to 42.2 A which does not 
seem to confirm the expected correlation between the EWs 
of 1 and Ha for Z CMa. It should be kept in mind though 
that the Ha emission peak contains an unknown component, 
which is emitted by the outflow region (Bailey 1998). The 
variation of the Ha emission EW is therefore partly due to 
variation in the outer atmosphere, whereas the emission of 
1 (8446.5 A) originates in the inner part of the atmosphere. 

Next to the Balmer emission lines, the Call (2) triplet 
lines have by far the strongest emission fluxes in the spec- 
tra of Z CMa. Although the three lines are blended by the 
Paschen lines P16, P15 and P14, the contributions from 
these lines to the flux are small (their EWs are about 5% of 
the total EW). The lines have P-Cygni profiles with outflow 
velocities around 500 km s -1 and half- widths of 107, 110 
and 132 km s _1 for 8498, 8542 and 8662 A respectively. The 
line shapes are very asymmetrical: the red wings extend to 
more than 8 A from the line center, which suggests that the 
lines are formed in a turbulent region. 

We confirm the presence of many emission lines from 
Fen of multiplets 27, 28, 32, 35, 36, 37, 38, 40, 42, 43, 46, 
48, 49, 55, 73, 74 and 199, several lines of Cm multiplet 30, 
43 and 44 and a few lines of Tin (20) and Sen (29). Apart 



from multiplet 73 in the red part of the spectrum, most of 
these lines are also given in the line list of the high state 
spectrum of Hessman et al. (1991). We discuss the variation 
of the Fen lines in Sect. 4.6. In the red part of the spectrum 
we also observed fairly strong emission of the lines of Fe I 
multiplet 60. 

4-2.2 Source of the broad absorption troughs of Hel and 
Ol 

The lines of Hel (7065, 6678 and 5876 A) (Fig. 5) and Ol 
(7773 and 6156 A) (Fig. 6) have relatively broad absorption 
features. They also show a weak emission at a small redshift 
(+65 km s _1 ) from the laboratory wavelength and a broad 
(—450 km s _1 ) blueshifted absorption trough and therefore 
have the shape of P Cygni profiles. There are several rea- 
sons to believe that these profiles are formed in the envelope 
of the primary Herbig BOe star: In Fig. 10 of the paper by 
Hessman et al. (1991) we see that in the post-eruption phase 
the 6678 A He 1 singlet line has disappeared and replaced by 
two Fe 1 lines with the double structure which is characteris- 
tic for the FUor component. The broad profile for the high 
state (in the same Fig. 10) shows only weak indications of 
these double lines, which suggests that during this state the 
continuum of the primary star is higher than that of the 
secondary FUor component (see also Sect. 3.2), so that the 
absorption line contribution from the primary (Be) compo- 
nent is dominating. In fact it was estimated by Whitney 
(1993) that in the high state the contribution of the FUor 
component to the visual continuum is only one third of the 
observed continuum and therefore one half of the contribu- 
tion by the primary Be star. 

The 7065 A and 5876 A Hel triplet lines in the Tauten- 
burg spectrum show the same profiles as that of the 6678 A 
Hel singlet line. Although the 6678 A line is outside the or- 
ders of the WHT spectrum, we can see from the comparison 
of the Tautenburg and WHT spectra that the profile of the 
5876 A line is the same in both spectra, so that the He 1 line 
did not change over a period of one month, in contrast to 
the neighbouring Nal profiles. Because of its high excitation 
energy (23 eV), the emission components in the Hel lines 
will originate from the hot neighbourhood of the stellar sur- 
face. The broad blue-shifted absorption components must be 
interpreted as a strong and fast outflow with a broad (per- 
haps turbulent) velocity spectrum. We also note here that 
no line of He 11 (e.g. 4686 A) was found in our optical spectra 
so that there are no direct indications for the presence of a 
chromosphere or transition region. 

The O 1 (7773 A) triplet has been observed in absorp- 
tion for several Herbig Ae/Be stars (Felenbok et al. 1988; 
Hamann & Persson 1992b). In Z CMa, this line shows indi- 
cations of emission components, redshifted by ~ 70 km s _1 
and the absorption components are very similar to that of 
the 6678 A He I line. This indicates that the line is not neces- 
sarily formed by the FUor component. It is also clear that in 
Jan. 2000 the continuum of the primary BOe star was much 
higher than that of the FUor at the time of the observation. 

Hamann & Persson (1992b) suggest that the Ol profile 
is formed in a highly turbulent envelope region. This may 
also be concluded from spectral surveys over many stars 
with different spectral and luminosity class (e.g. Slettebak 
1986; Jaschek et al. 1993). These authors show that the only 
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early B-type stars with Ol (7773 A) in absorption are the 
shell stars. In addition it has been shown by Faraggiana et 
al. (1988) that for this line EWs of the order of 2 A are only 
found for stars of spectral type earlier than Bl. For Z CMa 
the maximum velocity of the outflow is ~ 500 km s _1 in the 
Tautenburg spectrum. Because the excitation potential of 
the O I lines is only 10.7 eV, the inner border of the turbulent 
region were the emission of O I is formed will be further away 
from the stellar surface than that of He I (6678 A). The same 
profile as for the Ol (7773 A) triplet is found for the 6157 A 
triplet of O I, which, when in emission, cascades down by 
7773 A emission from the 3p 5 P state to the metastable 3s 5 S° 
state at Xi = 9.5 eV. 

4.2.3 The NalD, Call K and Kl lines 

The observed Na I D profiles (Fig. 5) consist of at least three 
components: a narrow absorption component, an emission 
component and a broad absorption component, composed 
of several blue-shifted lines. The narrow components of Di 
and D2 are shifted by ~ 27 km s _1 , which is close to the 
radial velocity of Z CMa. We assume that they are due to 
absorption by the foreground interstellar medium and the 
local CMa cloud. Their FWHM is ~ 24 km s _1 . The EWs 
of the narrow components (Table 4) can be used to estimate 
the Nal foreground column densities with the help of the 
Stromgren (1948) doublet method. By assuming a galactic 
abundance ratio 7V(H)/iV(Na) of 5 x 10 we can derive an 
average interstellar gas column and, with the empirical re- 
lation of Hobbs (1974), we can obtain a value for E(B — V) 
of 0.32, which is in line with the foreground excesses of the 
other stars in CMa Rl (papers I and II). 

In the high resolution spectra obtained at ESO (V ~ 
10. 2 m ) and Tautenburg (V ~ 9.1 m ), the emission compo- 
nents of the NalD lines are narrow and on both sides of 
the narrow absorption components. We suspect that they 
are the emission components (intersected by the interstellar 
absorption lines) of NalD P-Cygni profiles. 

The deepest and narrowest absorption component ex- 
tends from —30 to —200 km s _1 with respect to the in- 
terstellar component with a FWHM of 65 km s _1 . Accord- 
ing to Whitney et al., the deep absorption component of 
the D-lines is associated with the secondary (FUor) compo- 
nent. It causes a strong increase in polarization of the com- 
bined spectrum over this wavelength interval because the 
deep NalD absorption lines of the secondary source can- 
not dilute the polarized continuum emission of the primary. 
The deep absorption components around 5896 and 5886 A 
(F/F c w 0.2-0.3) will therefore almost certainly be formed 
around the secondary star. This is very similar to the Na I D 
line profiles of the FUors FU Ori and V1057 Cyg (see e.g. 
Welty et al. 1992, Figs. 6 and 7), where the maximum out- 
flow velocities of the absorption trough are —200 and —250 
km s _1 respectively. The (P Cyg) emission peaks are also 
usually seen for the Na I D profiles of FUors. At bright phases 
of Z CMa (e.g. the 'high state') the blue absorption wings 
(shortward 5886 A and 5890 A), however, may receive signif- 
icant contributions from the primary source, because then 
the continuum of the secondary star is lower than that of the 
primary. In Jan. and Feb. 2000 the absorption wings extend 
from —200 to —500 km s _1 . The fact that this new wing 
shows up simultaneously with the maximum in brightness 
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Figure 7. Fe II line profile variations in Z CMa between 1991 and 
2000. (a) Fen (46) 5991.38 A line profile variations, (b) the same 
for the Fen (74) 6238.38, 6247.56 A lines, (c) Fell (40) 6432.65 A 
line, (d) Fell 9997.58 line. Note the marked increase in Fen line 
strength in the Jan. and Feb. 2000 spectra. 

of the primary star suggests that this outflow component oc- 
curs in the outer atmosphere of the B0 III star. The strength 
of this wing varies within one month (Fig. 5). 

The K I (7698.96 A) line profile (not shown here) shows 
the opposite behaviour. In Jan. 2000 it consists of a (broad) 
emission component and a narrow one (FWHM ~ 12 
km s _1 ). If the absorption component is interstellar, its EW 
leads to a Kl abundance which is 56 times lower than the 
Nal abundance derived here. This is 2.8 times lower than 
the solar abundance of Kl. In the WHT spectra of Feb. 
2000 a second, broader (FWHM ~ 33 km s _1 ) absorption 
component has appeared, which is blue-shifted by about —33 
km s _1 with respect to the first component. This broad com- 
ponent, which could be formed in the outer disc or circum- 
stellar region is also present in the spectra of 1991 and 1996. 
It demonstrates again that during the outburst variability 
in the outer shell of the primary Be star occurs within a 
time-scale of less than one month. The Ca II K profile of the 
Feb. 2000 WHT spectrum (Fig. 8) shows an outflow veloc- 
ity pattern which is similar to that of the Na I D lines in the 
same spectrum, but completely different from those of the 
Fell (42) and the Balmer lines. This suggests that the CallK 
line is formed over the same outer envelope or disc region 
as NalD, which is in line with the comparable excitation 
potentials of the two lines. 



4.3 The spectra of Nov. /Dec. 1991 (V = 9.6 m ) 

4-3.1 The emission lines 

Although for Z CMa indications for spectral changes related 
with brightness variation have been reported previously in 
the literature (Covino et al. 1984; Finkenzeller & Jankovics 
1984; Hessman et al. 1991; Hamann & Persson 1992b; Welty 
et al. 1992; Teodorani et al. 1997; Garcia et al. 1999), a 
systematic comparison of red and blue spectra of Z CMa 
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Table 3. Lines detected in the blue spectra of Z CMa (bl = blend). Lines included in square brackets indicate forbidden transitions. 



A [A] 



Idcnt . 



WHT, 01/23/2000 



BTA, 11/22/1991 



WHT, 12/26/1996 







EW 


FWHM 


fout 


EW 


FWHM 






[A] 


[km s- 1 ] 


[km s- 1 ] 


[A] 


[km s -1 


3933.66 


Can (1) 


—0.96 


97 


-600 






4233.17 


Fc n (27) 


-0.44 


80 


-106 


-0.31 


130 


4278.13 


Fe II (32) 


—0.19 


88 




—0.17 


68 


4290.22 


Tin (41) 


-0.28 


53 


-107 


-0.08 


67 


4295.17 


Fei(41) + 


-0.28 


53 


-107 


-0.07 


49 


4295.90 


Tin (20) (bl) 












4306.67 


Fei(42) + 


-0.24 


44 


-98 


-0.08 


98 


4307.15 


Tin (41) (bl) 












4320.54 


Till (41) + 


-0.22 


67 


-79 


-0.07 


62 


4321.95 


Sen (15) (bl) 












4341.77 


H7 


-1.68 


158 


-704 


-0.48 


115 


4351.75 


Fc 11 (27) 


-0.56 


79 


-95 


-0.18 


115 


4358.36 


[Fen] 


-0.04 


53 




-0.03 




4368.25 


Fcl (41) 


-0.21 


61 


-61 


-0.19 


79 


4369.41 


Fen (28) 


—0.15 


61 




-0.09 


48 


4374.61 


Tin (93) + 


-0.26 


69 


-103 


-0.08 


73 


4375.69 


Sen (14) (bl) 












4384.32 


Fc 11 (32) 


-0.23 


64 


-95 






4385.38 


Fe 11 (27) 


-0.37 


60 




-0.16 


79 


4395.08 


Tin (19) 


-0.41 


74 


-145 


-0.16 


126 


4416.83 


Fc 11 (27) + 


-0.37 


60 




-0.17 


78 


4417.73 


Ti 11 (40) (bl) 












4468.49 


Tin (31) + 


-0.28 


51 


-114 


-0.14 


69 


4469.15 


Tin (18) (bl) 












4489.18 


Fc 11 (37) 


-0.49 


80 


-114 


-0.18 


115 


4491.41 


Fc 11 (37) 


—0.52 


57 




—0.14 


138 


4501.37 


Tin (31) 


-0.39 


59 


-118 


-0.08 


70 


4508.28 


Fc 11 (38) 


-0.49 


39 


-80 


-0.27 


95 


4515.34 


Fc 11 (37) 


—0.66 


69 


-80 


-0.18 


80 


4520.22 


Fc 11 (37) 


-0.39 


51 


-93 


-0.29 


100 


4522.63 


Fc 11 (38) 


-0.43 


63 




-0.4 


100 


4533.97 


Tin (50) + 


-0.34 


59 


-92 






4534.17 


Fc 11 (37) (bl) 












4541.52 


Fc 11 (38) 


-0.27 


59 








4549.47 


Fen (38) + 


-0.45 


56 




-0.32 


80 


4549.62 


Tin (82) (bl) 












4555.89 


Fen (37) 


-0.41 


61 




-0.21 


76 


4558.66 


Cm (44) 


-0.39 


56 




-0.26 


76 


4571.09 


Mgi(l) + 


—0.26 


51 


-81 


—0.21 


118 


4571.97 


Tin (82) (bl) 












4576.33 


Fc 11 (38) 


-0.47 


72 




-0.31 


134 


4583.83 


Fc 11 (38) 


-0.61 


77 


-122 


-0.51 


84 


4588.32 


Cm (44) 


-0.26 


55 




-0.19 


92 


4629.34 


Fc 11 (37) 


-0.62 


66 




-0.29 


112 


4634.11 


Cm (44) + 


-0.14 


44 




-0.08 


95 


4734.12 


Fc 11 (25) (bl) 












4708.66 


Tin (49) 


-0.08 


54 




-0.05 


61 


4731.44 


Fc 11 (43) 


-0.31 


59 




-0.11 


109 


4805.11 


Tin (92) 


-0.32 


52 


-87 


-0.04 


40 


4812.35 


Cr II (30) 


-0.06 


51 




-0.08 


64 


4814.53 


[Fen] 


-0.06 


55 




-0.08 


64 


4824.13 


Cm (30) 


-0.15 


55 




-0.04 


61 


4861.33 


H/3 


-5.92 


110 


-965 


-3.03 


123 


4876.41 


Cm (30) 


-0.11 


55 




-0.05 


47 


4923.92 


Fc 11 (42) 


-1.04 


85 


-186 


-0.52 


102 


4993.36 


Fc 11 (36) 


-0.39 


42 




-0.18 


78 


5018.43 


Fc 11 (42) 


-1.73 


123 


-230 


-0.95 


120 


5129.14 


Tin (86) 


-0.04 


33 




-0.04 


45 


5132.67 


Fc 11 (35) 


-0.33 


52 




-0.08 


67 


5158.78 


[Fen] 


-0.18 


82 




-0.13 


50 


5169.03 


Fc II (42) 


-1.91 


119 




-0.83 


77 



fout 

[km s- 1 ] 



EW 

[A] 



FWHM 

[km s- 1 ] 



fout 

[km s- 1 ] 



-1.64t 



71 



-650 



-760 



-54 



-48 
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Table 3. (Continued) 



A [A] Idcnt. WHT, 01/23/2000 BTA, 11/22/1991 WHT, 12/26/1996 



EW FWHM tw EW FWHM v out EW FWHM ?w 

[A] [kms- 1 ] [kms- 1 ] [A] [kms- 1 ] [km s" 1 ] [A] [km s^ 1 ] [km s" 1 ] 



5188.71 


1 1 II (70) 


—0.26 


44 




—0.11 


78 








5197.57 


TV t i / A f\ \ 

bell (49) 


—0.69 


66 




—0.51 


130 








5211.54 


1 1 II (103) 


—0.17 


73 




—0.11 


59 








5226.53 


11 II (7UJ 


—0.4/ 


66 




—0.06 


52 




—0.16 


83 


5234.62 


TV T T / A (\\ 

be II (49) 


—0.44 


80 




—0.23 


74 




—0.21 


59 


5237.37 


Cm (43) 


—0.19 


58 




—0.09 


81 




—0.09 


63 


5261.62 


[Fc n] 


—0.09 


58 




—0.05 


58 








5269.54 


Fe I (15) 


—0.45 


91 










—0.14 


60 


5275.99 


TV tt / A t~\ \ 

bell (49) 


—0.56 


58 




—0.25 


88 




—0.15 


65 


5279.88 


Cr II (43) 


—0.11 


65 




—0.04 


58 








5284.09 


TV „ / li \ 

le II (41) 


—0.81 


65 




—0.16 


73 




—0.07 


40 


5305.85 


Cr II (24) 


—0.14 


80 




—0.08 


40 








5316.61 


II (49 j + 


— 1.24 


72 




—0.64 


94 




—0.52 


91 


5316.78 


bell (48) (bl) 


















5325.56 


Tn„ tt / /i n\ 

II (49 j 


—0.26 


72 




—0.11 


58 




0.04 


36 


5328.04 


Fe I (15) 


—0.33 


67 










0.08 


45 


5333.65 


[Fc n] 


—0.07 


43 










—0.08 


60 


5336.81 


rp- tt / r-r\\ 

1 i ii (69) 


—0.35 


69 










0.11 


36 


5362.86 


TV t i / ylo\ 

1c II (48) 


—0.71 


33 




—0.37 


115 




—0.13 


43 


5371.51 


Fei (15) 


—0.17 


57 










—0.07 


50 


5397.13 


Fei (15) 


—0.21 


59 










—0.11 


47 


5405.78 


Fe I (15) 


—0.15 


56 










—0.06 


38 


5429.67 


Fe I (15) 


—0.18 


47 










—0.06 


46 


5432.98 


Fe ii (55) 


—0.17 


66 










—0.06 


53 


5434.53 


Fei (15) 


-0.08 


60 










-0.02 


18 


5446.92 


Fei (15) 


-0.11 


47 










-0.05 


28 


5455.62 


Fei (15) 


-0.14 


51 










-0.04 


42 


5497.53 


Fei (15) 


-0.11 


65 










-0.07 


46 


5501.48 


Fei (15) 


-0.06 


46 










-0.03 


35 


5506.79 


Fei (15) 


-0.06 


69 










-0.04 


46 


5534.86 


Fc ii (55) 


-0.53 


64 










-0.11 


48 


5577.34 


[Oi] 


-0.05 


73 










-0.15 


37 


5591.37 


Fc ii (55) 


-0.74 


68 










-0.09 


36 


5754.64 


[Nil] 


-0.05 


30 










-0.02 


40 


6562.81 


Ho 


-44.1 


191 


-920 


-16.8 


200 


-850: 


-18.6+ 


96 



t Observed with ESO CAT, 12/16/1996. 

in its "low" and "high"-state has not been possible so far. 
Using previously unpublished spectra from Z CMa in its 
"low state", taken from the ESO data-archive, as well as 
our new data of Z CMa after its Jan. 2000 outburst, we can 
now compare spectra of Z CMa in its two states over the 
full spectral range between 4200 and 10000 A. 

From this comparison we draw the following conclu- 
sions: (1) The EWs of the Paschen and Call (2) emission 
components in the Jan. 2000 spectra are a factor of 2.0-2.5 
larger than the corresponding EWs in the Dec. 1991 spectra. 
We discuss the consequences of this in Sect. 4.5. The half- 
widths of the Call lines are comparable to those of Jan. 2000 
but the outflow velocities are 140-180 km s _1 lower than in 
2000. (2) The EW of the Oi (8446 A) emission line does not 
seem to increase significantly between the observations in 
1991 and 2000 but it is ~ 16% larger in Dec. 1987 (Hamann 
& Persson 1992a) than in Jan. 2000. If this line is excited 
by Ly/3 pumping, it indicates that it is formed at a location 
where the Ly/3 fluxes are not significantly different for Dec. 
1991 and Jan. 2000. (3) The weak Nl emission lines at 8728 



and 8683 A, seen in the Jan. 2000 spectra, are also present 
in the 1991 NTT spectrum. (4) All obtained spectra contain 
many emission lines of Fell, Tin and Cm. The analysis of 
the Fell lines is given in Sect. 4.6, together with that of 
the outburst spectra of 1987 and 2000. We also observe the 
emission lines of Fei (60). (5) In contrast to the other Fell 
emission lines, the Fen (9997 A) line, which is extremely 
strong in the spectra of Jan. /Feb. 2000, is not present in the 
NTT spectrum of 21 Dec. 1991. This line is the strongest line 
in the near-IR multiplet b 4 G — z 4 F° from which four lines 
are observed in the spectra of Jan. /Feb. 2000 (Table 4) and 
in which the 9997 A line is by far the strongest transition. 
The lines may be pumped directly from Lya (fluorescence), 
by cascading from other Lya-pumped levels, or collisionally 
pumped from the metastable a 4 G level. These three possi- 
ble mechanisms correspond with different excitations in the 
UV (Rodriguez- Ardila et al 2002). Without high resolution 
UV spectra of Z CMa the method of formation cannot be 
selected. 
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Table 4. Lines detected in the red spectra of Z CMa (bl = blend, s = single, tr = triplet). Lines included in square brackets are forbidden 
transitions. 



\ r X 1 
A [A] 


Ident. 


Tautenburg, 01/23/2000 


WHT, 02/18/2000 


NTT, 12/21/1991 


WHT, 12/26/1996 






EW 


FWHM 


fout 


EW 


FWHM v out 


EW 


FWHM 




EW 


FWHM j; out 






[A] 


[km s- 1 ] 


[km s- 1 ] 


[A] 


[km s- 1 ] [km s- 1 ] 


[A] 


[km s- 1 ] 


[km s- 1 ] 


[A] 


[kms- 1 ] [kms- 1 ] 


5657.87 


Sen (29) 


— U.Zo 


QQ 
OO 














A AQ 

U.Uo 


A Q 

4o 


5875.62 


Hel(tr) 


O OQ 

Z.Zo 


1 AQ 

iyo 


A OA 

— 4zU 


A Q/1 

U.o4 


1 OQ Q77 

iyo — o * / 












5889.95 


Nal(Di) + 


a oo 


1 OQ 

IZo 


ftftc; 
— 000 


A A1 

4.U1 


1 QO ft C^A 

ioy —oou 


O AQ 

Z.Uo 


1 OQ 

Izo 


c;i 7 
— 01 ( 


Q A 1 

o.41 


1 qa cot; 
loU — OZo 


5895.92 


Nal(D 2 )(bl) 






















5991.38 


Fen (46) 


a 7ft 
— U. ( 


11 




A fi7 
— U.O ( 


7/1 
( 4 


A O^ 


fi/1 
04 




A AO 

—u.uy 


/1Q 
40 


6084.11 


Fen (46) 


a oo 
— U.zy 


1 AO 

iuy 




A Q7 
— U.O / 


1 AA 
1UU 


A A7 
— U.U ( 


ftQ 

Oo 








6113.33 


Fen (46) 


a Aft 

— U.UO 


tr 
00 








A A/1 
— U.U4 


R7 
( 








6116.61 


Fen (46) 


a n7 
— U.U i 


00 


















6129.73 


Fen (46) 


— U.Uo 


ftft 

OO 














O AO 

— u.uz 


QQ 
OO 


6136.33 


Fen (169) 


A 1 Q 

— U.lo 


A O 

4z 


















6147.74 


Fen (74) + 


— l.Zo 






A QQ 
— U.OO 




A OQ 
— U.ZO 






A Ol 

— U.Zl 




6149.25 


Fen (74) (bl) 






















6156.41 


Oi(10) 


a oo 
u.zy 






A OO 

u.zy 


OQA 
ZoU 












6191.75 


Fci(169) 


A Q1 

— U.ol 


ftQ 

Oo 




A QQ 

— U.oo 


ft1 
01 


A AO 

— u.uy 


f;A 

OU 




A AQ 

— U.Uo 


Q er 
oO 


6238.39 


Fe II (74) 


A OO 

— u.yy 


O/l 

y4 


QftA 
— oOU 


A ftf; 

— U.OO 


Qft 
OO 


— U.lO 


7Q 
/ o 




A AO 
— U.UZ 


Q1 
ol 


6247.56 


Fe II (74) 


i no 


on 

oy 


Q/1 A 

— o4U 


A Q/1 

— U.o4 


7Q 
( O 


A Q A 

— U.o4 


77 








6252.56 


Fei(169) 


A Q1 

— U.ol 


Oo 




A 1 ^ 
— U.lO 


Oo 








A AQ 
— U.UO 


7A 
I U 


6300.31 


[Oi](l) 


A AA 
— U.44 






A A fit 
— 4.40 1 




A Qft 

—U.oO 






O AQ 
— Z.OO 




6318.13 


Fei(168) 


A 07 
— U.Z ( 


A A 

4U 




A 1 Q 

—U.lo 


It A 

04 


A Aft 

—U.UO 


OZ 








6331.95 


Fen (199) 


—0.08 


88 




—0.06 


60 


—0.08 


a a 
04 








6347.09 


Si ii (2) 


A 1 1 

— U.ll 




















6363.79 


[Oi](l) 


A fl7 
— U.U ( 






A 1 1 

— U.ll 


00 


A A7 
— U.U / 


a a 
04 




A 1 O 
— U.IZ 


00 


6369.46 


Fe ii (40) 


A 1 Q 

— U.lo 


fti 

01 




— U.lO 


A O 

4y 


Alt; 
— U. 10 


c;o 
OZ 








6371.36 


Si n (2) 


— U.lO 


A O 

4z 


















6393.61 


Fei(168) 


A Ol 

— U.zl 


Q/1 

o4 




A 1 Q 

— U.lo 


t;ft 
00 


A Aft 

—U.UO 


c;o 
OZ 




A AQ 

— U.Uo 


1 ft 
10 


6407.29 


Fen (74) 


A 1 Q 

— U.lo 


1 AO 

iuy 




A 1 O 

— U.IZ 


1 A/1 

1U4 


A A7 
— U.U ( 


KA 

OU 








6416.91 


Fen (74) 


-0.69 


85 




-0.37 


69 


-0.2 


61 




-0.09 




6432.65 


Fen (40) 


1 QO 
— l.OZ 


7fi 
(0 




1 Oft 

— 1 .zo 


1 O 


A 07 
— U.Z ( 


00 




A 1 7 
— U. 1 ( 


A1 
4 ( 


6456.38 


Fen (74) 


1 Qft 
— l.OO 


Q K 

oo 








A Qt; 

—U.oO 


/ 








6482.15 


Fen (199) 


A 1 7 
— U.l / 


7ft 

to 




A Oft 

— u.zu 


51 












6516.81 


Fe II (40) 


1 7Q 
— 1. ( O 


Q A 

o4 


QQA 

— ooU 


1 OQ 

— l.Zo 


ftA 

OU 


A C^Q 

— U.Oo 


7/1 
/ 4 




A OA 

— u.zy 


A Q 
4o 


6548.04 


[Nn] 


—0.03 


"1 7 
I ( 


















6562.81 


Ha 


AO A 


OAA 
ZUU 


1 1 AA 
— 11UU 


A A 1 
— 44. 1 


1 Ol OOA 

iyi — yzu 


1 1 o 

— 1 1 .y 


1 OA 

IZU 


1 A/1 A 

— 1U4U 


OQ A 

— zo.4 


Oft 1 AAA 

yo — iuuu 


6583.46 


[Nn] 


a no 
— U.UZ 


1 Q 

lo 




A AO 

— U.Uz 


oo 

zz 


A AO 

—U.UZ 


OQ 

Zo 




A AQ 

— U.Uo 


1 7 
1 ( 


6678.15 


Hel(s) 
Hcl(s) 


A fi7 
U.O / 

a riQ 
— U.Uo 


QAA 

oUU 








A Aft 
U.UO 


c;q 
Oo 








6707.81 


Lii(l) 


A Aft 
U.UD 


OC 

Zo 








A ac; 
U.UO 


OQ 

Zo 




A AA 

u.uy 




6717.24 


[S n] 


A 1 Q 

— U.lo 


OA 

zu 




n 1 fi 
— U. 10 


ftA 

OU 








A AQ 

— U.Uo 




6726.28 


Oi(2) 


-0.02 


75 


-200 
















6731.31 


[S ii] 


A A7 

— U.U ( 


OA 

ZU 




A 1 O 

— U.IZ 




A AQ 

—U.Uo 






a At; 
— U.UO 




6757.21 


Si (8) 


A A1 
— U.U1 


QQ 
OO 


















7005.89 


Si I (60) 


—0.05 


33 


— 110 


—0.04 


58 












7017.67 


Si I (51) 


A A/1 

— U.U4 


A A 
44 


1 1 A 

— 11U 
















7065.19 


He I (tr) 


A A A 
— U.44 


OQ A 

ZoU 


A QA 
— 4oU 


A A7 
— U.U ( 


con 

oyu 


n A/i 
<U.U4 










7155.16 


[Fen] 








A 1 A 
— U.1U 


f^A 

ou 


A 1 A 

— U. 1U 


OO 








7244.85 


Si I (51) 


A A/1 

— U.U4 


A A 
44 


1 1 A 
— 11U 






A 1 A 

— U. 1U 










7281.35 


Hci(s) 
Hci(s) 


A ftf; 
— U.OO 

A 1 A 

U.14 


QA 
OU 








A A7 
— U.U ( 










7291.46 


[Can] (1) 


A £1 

— U.ol 


c;ft 
00 


77 






A 01 
— U.Zl 


fiQ 

Oo 




A Ol 

— U.Zl 


Ol 

Zl 


7310.21 


Fen (73) 


A Qf^ 
— U.oO 


<.yo 








A 1 1 
— U. 11 


fiO 

oy 








7320.69 


Fen (73) 


-0.27 


55 








-0.11 


66 








7323.88 


[Can] (1) 


— U.40 


oo 


11 
— ( ( 






fl 11 
— U.ll 


OO 




A 1 A 

— U. 14 


Q A 

oU 


7377.91 


[Mil] (2) 


-0.04 


24 


-60 


-0.01 


27 












7411.87 


[Nin](2) 


-0.08 


64 








-0.03 


31 




-0.05 


55 


7415.96 


S I (22) 


-0.19 


100 


-100 


-0.15 


87 


-0.04 


51 








7449.34 


Fen (73) 


-0.32 


86 




-0.28 


78 


-0.15 


90 








7462.34 


Fen (73) 


-0.71 


68 




-0.57 


72 


-0.14 


48 








7515.84 


Fen (73) 


-0.25 


60 




-0.22 


75 


-0.08 


47 




-0.03 


71 


7698.96 


Kl(l) 
Kl(l) 


0.08 
-0.08 


12 




0.32 
-0.09 


15 


0.20 
-0.05 






0.17 
-0.08 




7711.73 


Fen (73) 


-1.11 


82 




-0.9 


76 


-0.25 


66 




-0.06 


50 
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Table 4. (Continued) 



A [A] Ident. Tautcnburg, 01/23/2000 WHT, 02/18/2000 NTT, 12/21/1991 WHT, 12/26/1996 

EW FWHM v out EW FWHM v out EW FWHM « out EW FWHM v ou t 

[A] [km s- 1 ] [km s- 1 ] [A] [km s" 1 ] [km s" 1 ] [A] [km s" 1 ] [km s" 1 ] [A] [km s" 1 ] [km s" 1 ] 



7774.17 


Ol(l)(tr) 
Oi(l) 


2.79 
—0.08 


340 


450 


3.12 
—0.19 


300 


—250 


0.11 










7877.51 


TV T TT /o\ 

Mg 11 (8) 


—0.18 


80 


—304 


















701^/* no 

7896.38 


TV if T T / O \ 

Mgii(8) 


—0.23 




—304 








—0.03 










7918.38 


Si 1 (57) 


—0.08 


66 




















7932.35 


Si 1 (57) 


—0.08 


66 




















■70/1 A 1 1 

/ y44. 11 


Ol I (0 / ) 


— U.U / 


00 










— U.UO 










o nnn i > i 

8000.12 


rn TT i /-1 \ 

[Cm] (1) 


—0.21 


60 




















8U12.27 


S 1 (10) 


—0.04 


44 




















8047.81 


Mg 1 (0) 


—0.17 


66 


—250 














—0.05 


32 


8201.71 


Can (13) 


—0.15 


60 




—0.18 


87 




—0.05 






—0.05: 




oni o nn 

8213.99 


TV If T T / 7\ 

Mg II (7) 


—0.03 


60 




—0.07 


78 














8234. o4 


Mgn (7) 


—0.14 


60 




—0.22 


82 




—0.14 


42 








8248.81 


Can (13) 


—0.27 


64 


—550 


—0.27 


83 


—250 












O O O 7 "1 "1 

8327.11 


Fc 1 (60) 


—0.26 


60 




—0.27 


76 




—0.12 


51 








8345.45 


P23 + 


—0.53 


120 










—0.11 


39 








8346.13 


TV X f A f\\ /l 1 \ 

Mgi (40) (bl) 
























OO09.UU 


POO 1 


— U.4 / 


Izo 










— U. ly 


1 1 n 
11U 








8361.69 


He 1 (bl) 
























8374.48 


P21 


—0.44 


80 




















O Q O 7 70 
OOO /. (8 


1c I {bl)) 


—0.64 


80 










—0.19 


58 




—0.19 


60 


8392.40 


P20 


—0.63 


100 




















Olio in 

8413.32 


P19 


—0.54 


110 


—440 


















8446.51 


Oi(4) 


— 1.37 


107 




1.82 


165 


—270 


— 1.31 


120 








040 I .ZD 


pi *7{U1\ 


— U.Oo 


yo 


— 4oU 


— U.Oo 


111 
111 


— ZOO 


n 17 
— u. 1 1 


1 








84o8.42 


Fc 1 (60) 
























8498.02 


Can (2) 


— 12.6 


107 




— 10.86 


83 


—430 


—5.51 


93 


—250 


—4.76 


101 


8502.41 


P16(bl) 
























8514.10 


Fe 1 (60) 


—0.48 


62 




—0.38 


66 




—0.07 


60 




—0.05 


44 


8542.10 


Can (2) + 


— 14.9: 


110 


—500 


— 11.3 


81 


—460 


—5.82 


109 


—320 


—4.95 


79 


8545.38 


P15(bl) 
























orno on 

8598.39 


P14 


—0.61 


89 


—420 


—0.37 


83 




—0.19 


31 








8629.24 


N 1 (2) 


—0.17 






















8662.14 


Can (2) + 


— 11.8 


132 


—440 


— 14.9 


73 


—500 


—5.26 


106 


—277 


—3.59 


109 


8665.02 


P13(bl) 
























8683.41 


Ni(l) 


—0.06 


55 










—0.03 


10 








8688.64 


Fe 1 (60) 


—0.65 


76 










—0.12 


50 




—0.08 


50 


8728.89 


Ni(l) 


—0.16 


56 










—0.18 


20 








O 7 7* n A 7 

8750.47 


P12 


—0.64 


93 


—430 


—0.29 


87 




—0.32 


64 








8806.74 


Mgi(7) 


—0.84 


111 


—325 


—0.83 


101 




—0.06 


35 




—0.04 


20 


O O O A "1 "1 

8824.11 


Fe 1 (60) 


—0.51 


110 




—0.48 


66 




—0.08 


66 




—0.05 




OO/Tl 70 

8862.78 


Pll 


—0.64 


85 


—433 


—0.51 


85 




—0.28 


52 








8912.07 


Can 


—0.58 


106 


—415 








—0.05 


40 








o nn 7 o / ■ 

8927.36 


Can 


—0.82 


106 


—420 








—0.06 


51 








nn "i ■! n "i 

9014.91 


P10 














0.32 


54 








t \ i \ / ' - > n 7 

9063.27 


He 1 (77) 


—0.68 


150 


—435 


















nnnn nn 

9229.02 


P9 


—0.61 


95 


—410 








0.49 


66 








9545.97 


P8 


— U.oZ 


on 
yu 


a An 

— 44U 


















9854.74 


Can 


-0.43 


106 




















9997.57 


Fen 


-2.02 


110 


-360 


-1.71 


110 


-380 


0.00 






0.00 




10049.38 


P7 


-1.1 


88 


-420 


-0.85 


84 


-342 


0.00 






0.00 




10370.51 


[Sir] 








-0.08 
















10397.74 


[Ni] 








-0.16 
















10459.79 


[Nin] 








-0.31 
















10501.50 


Fen 








-0.09 
















10525.14 


Fen 








-0.16 
















10931.83 


Pi 








-0.23 
















10937.83 P6 (P7) 








-0.10 

















-282 



-307 



-430 



t Observed with ESQ 3.6 m, 03/31/2002. 
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4-3.2 The absorption lines 

The strong and broad He I (7065, 6678 and 5876 A) and 
Ol (7773 and 6156 A) absorption troughs in the spectra of 
Jan. /Feb. 2000 are much less prominent in the spectra of 
Dec. 1991. For these lines the same variation was observed 
by Hessman et al. (1991) between the 'high state' and the 
'post-eruption state'. It suggests that the contribution by 
the FUor dominates the local continuum in Dec. 1991. This 
is probably also the reason, why the He I (5876 A) absorption 
is not present in the NalD spectrum of Dec. 1991 and Dec. 
1996 (Fig. 5). We also find a corresponding reduction in the 
strength of the He I 7068 A line, together with the appear- 
ance of narrow absorption lines of Fe I at 7068 A and of Fe I, 
Til, Mm at 7069-7070 A and Fel at 7071 A. The shape of 
the O I (7773 A) profile is shallow around 7772 A and shows 
a very small emission peak around 7776 A. The blue wing 
of the 5889.95 NalDi line in the 2000, 1996 and 1991 spec- 
tra changes in shape. The first absorption component (v < 
130 km s _1 ) increases in depth between 1991 and Feb. 2000, 
but decreases between 1996 and 1991. However, the Jan. 
2000 spectrum shows a second absorption wing with v < 
380 km s~ , which suggests that both lines are formed in 
the variable outer parts of the envelope or disc atmosphere. 
The Kl (7699 A) line profile is similar to that of the WHT 
spectra of Feb. 2000. 

4.4 The spectra of Dec. 1996 (V = 10. 2 m ) 

4- 4-1 The emission lines 

The red WHT spectrum of Dec. 1996, the ESO echelle pro- 
files and the Ha profile of Bailey (1998) have been obtained 
during the deep photometric minimum (V ~ 10. 2 m ) between 
1996 and 1999 (Fig. 1). Parts of similar spectra, taken in 
Jan. 1996 and Jan. 1997 have been published by Teodorani 
ct al. (1997) and Chochol et al. (1998) respectively. Our 
low resolution spectra of Dec. 1992 also correspond to this 
brightness. Since the brightness is changing slowly during 
these years, we have no direct indications for the relative 
spectral contributions of the components of Z CMa for this 
phase. However, some hints can be given by comparing these 
spectra with those of 1991, 2000 and Feb. 1987 (see next sec- 
tion). (1) The ESO profile of Ha of Dec. 12 appears similar 
to the WHT profiles of Dec. 26. The EW of the Ha emission 
has increased slightly with respect to that of 1991. However, 
Chochol et al. (1998) have found that the daily variations of 
the Ha emission can be considerable so that small changes 
over long periods do not permit us to draw any conclusions. 
We also note here that the Paschen emission lines are not 
detected in the 1996 spectrum, so that the Call (2) lines do 
not have to be corrected for the contribution by the Paschen 
lines. (2) The ESO profiles of the Can (2) lines at 8498 and 
8542 A are very similar to those of the WHT spectrum of 
Dec. 26. The EWs of the strong emission lines of Call (2) 
and Fell are lower by a factor of ~ 0.6 than the correspond- 
ing EWs measured from the 1991 spectrum. We discuss the 
consequences for the radii of the emission region in sections 

4.5 and 4.6. The profiles of the Can lines are similar to those 
of Dec. 1991. (3) The EW of the Fel (60) line at 8387 A is 
reduced by a factor of 0.85. (4) The EWs of the forbid- 
den emission line of O I (6300 A) have gradually increased 
since 1987. For [Fen] 7155 A and the [Call] lines at 7291 



and 7324 A there also may be an increase with respect to 
the 1991 data. A probable explanation for this rise is given 
in Sect. 4.6. (5) The EW of the CanK emission varies less 
than those of the Call (2) triplet. It even slightly decreases 
between 1996 and the outburst stages of 1987 and 2000. This 
can possibly be due to the fact that the Call K line is formed 
in the outer envelope, in contrast to the Call (2) disc lines. 

4-4-2 The absorption lines 

(1) Although the He I (7065 A) line has not been observed 
in the high state of Feb. 1987 and although we observe only 
the blue half of the He I (6678 A) line in the 1996 spectrum, 
we can compare the remaining profiles of these lines in the 
available spectra. We find that the EWs of these lines are 
decreasing with decreasing brightness, which indicates a de- 
creasing contribution of the primary BOe star. We can ex- 
trapolate this increasing trend and find that in Dec. 1996 the 
contribution of the primary is at most 10%. We would ex- 
pect a corresponding increase of the EW for the Li I (6707 A) 
line (which is only contributed by the FUor component), but 
the relatively large errors in line flux due to the weakness 
of this line in the four spectra prohibit drawing any reliable 
conclusions regarding an increase in its EW. The decreasing 
trend is also found for the He I (5876 A) absorption line. In 
the 1996 spectrum this line is too weak to be distinguish- 
able in the neighbourhood of the NalD absorption doublet. 

(2) The Ol (7773 A) has not been observed in this spec- 
trum (it lies in the gap between order 36 and 37). (3) The 
NalD profiles of the ESO and WHT spectra of Dec. 1996 
are similar, but the WHT spectrum of Dec. 26 has a higher 
resolution and shows that the blue components of the lines 
are double with mean velocities of —53 and —128 km s _1 . 
The NalD and Kl profiles do not differ much from those 
observed in 1991. The comparison of our 1996 CanK pro- 
file can be made only with those of the outburst spectra 
of 1987 and 2000. Our 1991 profile has a higher emission 
component than the other two profiles. The blue absorption 
wings of the three profiles are very similar, but in the 2000 
spectrum we observe two narrow absorption components at 
—305 and —460 km s~ x , which are difficult to find in the 
other two spectra, because of the noise. 

4.5 The Call (2) emission lines 

In the four high resolution spectra described above, the 
Ca n (2) triplet emission lines at 8498, 8542 and 8662 A are 
(apart from Ha) the strongest lines. Following Hamann & 
Persson (1992a, b) we shall use our data on these lines as 
a diagnostic tool to obtain some information about the ex- 
tended envelope around Z CMa. The profiles of this triplet 
in the spectra of Feb. 2000, Dec. 1991 and Dec. 1996 are 
given in Fig. 8. It shows that they have P-Cygni shapes, 
as is characteristic for many Herbig and T Tauri stars with 
outflow. As noted by Hamann & Persson (1992a, b) other 
characteristics of this group are: (a) a tendency for weaker 
emission lines to be narrower (in order Ha, Call (2), Ol, 
Fell), (b) the Can (2) peak- flux decreasing in order of in- 
creasing wavelength and the EW of the 8542 Ca II line being 
the largest of the triplet EWs and (c) the presence of Fe I 
emission lines. 
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Figure 8. Call line profile variations in Z CMa between 1991 and 2000. From left to right: (a) the CanK line at 3933.66 A (this spectral 
range includes Fell 3938.29 A, seen in emission in the Feb. 2000 spectrum), (b) the 8498.03 and 8542.09 A line of the red Call (2) triplet 
(this spectral range also includes the 8514.08 A Fel line, and the 8467.8 A P17/Fcl blend), and (c) the 8662. 14 A line of the Can (2) 
triplet (blended with P13). Again note the marked increase in line strength in the Jan. and Feb. 2000 spectra. 



Table 5. Can (2) line-ratios relative to the 8542 A line in the 
spectra of Z CMa, V645 Cyg, FU Ori and V1057 Cyg. 



EW Peak Height 



8498/8542 8662/8542 8498/8542 8662/8542 



Z CMa (Jan. 2000) 


0.85 


0.79 


1.44 


0.73 


Z CMa (Feb. 2000) 


0.96 




1.06 




Z CMa (Dec. 1991) 


0.95 


0.90 


1.18 


0.94 


Z CMa (Dec. 1996) 


0.96 


0.73 


1.10 


0.96 


Z CMa (Nov. 1988) + 


0.87 


0.79 


1.00 


0.86 


V645 Cyg* 


0.86 


0.95 


1.18 


0.91 


FU Orit 


0.59 


0.88 


0.67 


0.89 


V1057 Cygt 


0.65 


0.93 


0.89 


0.92 



^Determined from Welty et al. (1992). ■^Determined from 
Hamann & Pcrsson (1989). 



In Table 5 we give the peak and EW ratios of the triplet 
lines for the five spectra of Z CMa and for those of FU Ori 
and V1057 Cyg. It is clear that the ratios of both FUors 
are different from those of Z CMa and V645 Cyg. Our main 
argument that the observed Ca II (2) emission is produced 
by the B star component of Z CMa is that the EWs of the 
Ca II (2) emission lines strongly increase with the increas- 
ing continuum contribution by the B-star (10%, 25%, 59%). 
However, in the spectrum of 1996, where the contribution of 
the BO-star is only about 10%, Fig. 8 shows a component in 
the red wing of the Call emission profiles, which could be 
due to the FUor component. We also find a difference in the 
EWs of the 8540 A and 8498 A Call (2) emission lines (both 
—3.08 A) from the ESO CAT observations, obtained on Dec. 
15, 1996, and those (-4.95 and -4.76 A) observed with the 
WHT on Dec. 26, 1996. This is consistent with a smaller 
FWHM (93 km s" 1 for both lines) on Dec. 15 than those 
from Dec. 26 (97 and 101 km s™ 1 ). These differences could 
be caused by a variation in the disc of the B-star and/or in 
the FUor disc. 

In our spectra of 2000 and 1991 the (extinction cor- 
rected) ratio of the total flux in the Ca II doublet at 8912.7 A 
and 8927.4 A to the total flux of the Call (2) triplet is 0.036. 
Ferland & Persson (1989) have estimated that in the tem- 
perature range 5,000-10,000 K, (where Can is the dominat- 
ing ionization stage) it then follows that n e < 10 11 cm" 3 . 



Because the EW ratio of the triplet members is quite dif- 
ferent from the ratio 0.11:1.0:0.56 for the values of \og(gf) 
(which should be the ratio if the lines were optically thin) 
it has been recognised (Polidan & Peters 1976) that the 
Ca II triplet emission lines are optically thick and saturated. 
Hamann & Persson (1992b) showed that for t = 1 a min- 
imum value for the Can column density can be given by 
iV(Call) = 2.5 x 10 15 x {v tulb /(50 km s" 1 )) cm" 2 pro- 
vided that the line-width is determined by turbulence. For 
Z CMa this means that A r (Call) > 2.5 x 10 15 cm" 2 and 
therefore N(E) > 1 x 10 21 cm -2 . The condition for satura- 
tion is r x n £ > 10 13 cm" 3 , so that the optical thickness of 
the line formation region should be r ^ 100 and therefore 
AT(Can) and A?"(H) much larger. 

Hamann & Persson (1992b) estimate the size of the Call 
emitting region (Rcmi/ R*) from the value of the surface flux 
of the 8542 A Call line, defined by (L nne / X*)crT c 4 ff . If this 
line surface flux just equals ttxAvxBv(T cx ), where Av is the 
full width at half max. of the line and T ex is the excitation 
temperature, they assume that the line emitting region just 
covers the stellar surface, while larger surface fluxes corre- 
spond with in proportion larger surfaces. The authors have 
calculated the latter function for T ox = 5000 K, v — 25 and 
250 km s" 1 and T ox = 10,000 K and v = 250 km s" 1 and 
compared (their Fig. 9) the line surface fluxes, derived from 
the observations of the 8542 A line for various Herbig and T 
Tauri stars, with the three calculated limits. For their data 
of Z CMa (26-28 Dec. 1987) they assumed that it is a single 
star with spectral type F5. Then the derived 8542 A line sur- 
face flux is equal to the T ex = 5000 K, v = 250 km s" 1 limit, 
so that the Ca II emitting region just covers a star with this 

c X • 

From our new data we derived that the Call emission 
is formed in the neighbourhood of the B0 star component 
of Z CMa. We can scale the line surface flux from that of 
Hamann & Persson (1992b) by noting that in Jan. 2000 our 
estimate of the extinction corrected (R cs = 4.2) continuum 
flux at 8542 A for the B star component is 0.7 times the 
corresponding flux of the single F5 star in Dec. 1987. On 
the other hand the EW of the 8542 A line in Jan. 2000 is 
3.4 times its value, measured in Dec. 1987 and the T e g of 
a B0III star (~ 30,000 K) is a factor 4.4 larger than that 
of an F5I (~ 6900 K) star. In this way we obtain for Jan. 
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2000 a line surface flux of 1.71 x 10 10 erg cm' 2 s _1 . This 
implies a radius of ~ 29.2 R*. Fig. 10 of Hamann & Pers- 
son (1992b) then shows that this radius corresponds to the 
Stromgren radius for a BO star if n e ~ 8 x 10 10 cm -3 , 
which is consistent with the value of n e estimated from the 
Ca II doublet and triplet observations above. Ferland & Rees 
(1988) and Ferland & Persson (1989) have made calculations 
of the photoionization in the circumstellar region and have 
shown that at high n e the Call emission can be produced 
just beyond the Stromgren radius. Although in this region 
the dominating ion is Cam, the emissivity of Call appears 
to be high there because of the high temperatures. For early 
type B stars the emission of Call can arise from extended 
high density regions from close to the photosphere up to the 
Stromgren radius, e.g. in and around circumstellar discs. 

We can obtain an estimate of the radius of the Ca II disc 
in Nov./Dec. 1991 by noting that for this epoch the EW 
of the 8542 A line is 2.2 times smaller and the extinction 
corrected (R cs = 4.2) continuum flux near the line is 3.4 
times smaller than the corresponding values for Jan. /Feb. 
2000. This means that the line surface flux in Nov./Dec. 
1991 will be 0.13 times the line surface flux in Jan. /Feb. 
2000, so that (if T ex is the same) it becomes 2.3 x 10 10 erg 
cm -2 s _1 , which implies an outer disc radius of ~ 10.7 R* 
in Nov./Dec. 1991. In the same way we can use the data of 
Dec. 1996. For this epoch we estimated that the observed 
continuum flux of the B0 component is about 10% of the 
total observed flux. From this we derive that the extinction 
corrected (R v = 4.2) flux near the 8542 A line is 4.17 x 
10 -12 erg cm -2 s , which is 2.3 times smaller than the 
corresponding flux in Jan. /Feb. 2000. The EW of the 8542 A 
emission in the Dec. 1996 spectrum is 3.7 times smaller than 
the corresponding EW in the Jan/Feb. spectra, so that the 
line surface flux in Dec. 1996 ~ 2.0 x 10 9 erg cr 
(if T cx is the same) the outer disc radius will be 
Dec. 1996. We conclude that if T ox did not change, the outer 
radius of the Ca II (2) emission region increased during the 
'outburst' of Jan 2000 with a factor of roughly 3. If T ox rises 
the increase in the outer radius will be smaller. 



4.6 The Fen emission line spectra 

4-6.1 Permitted Fell lines 

The EWs of the Fell emission lines in the spectra of Feb. 
1987, Nov/Dec. 1991 and Jan/Feb. 2000 have been mea- 
sured and are listed in Table 3 and 4. We analysed these 
line spectra in a statistical way: the self-absorption curve 
(SAC) method (see Appendix), which is a Curve of Growth 
version for emission lines. Strictly speaking we should cor- 
rect the emission EWs for the line absorption, which under- 
lies each emission profile. This can be done for a number of 
lines, which are observed in absorption in spectra of low ac- 
tivity phases (Finkenzeller & Jankovics 1984; Covino et al. 
1984; Welty et al. 1992). However, for the majority of the 
lines these absorption components are weak and difficult to 
estimate accurately. In addition part of the absorption com- 
ponents at low brightness phases may be contributed by the 
FUor component in Z CMa. Therefore we have neglected 
this correction for all lines. 

Because the emission lines contribute very little (at 
most a few per cent) to the continuum, we can use the 
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Figure 9. Self-absorption curves (SACs) for Fen lines in the 
Z CMa spectra from Feb. 1987 (top), Feb. 2000 (middle), and 
Dec. 1991 (bottom). The + symbols indicate the optically thick 
Fen lines, the © symbols indicate the group of optically thin 
Fen lines, and the thick drawn lines are linear least square fits 
to each of these groups. The thin drawn line is a linear least 
square fit to both Feu line groups. Typical observational errors 
are comparable to the size of the plot symbols. The thick crosses 
(+) in the upper left corners of the middle (SAC 2000) and lower 
panel (SAC 1991) mark the positions of the forbidden [Fell] lines. 
They are not included in the least square fits. The caption of each 
panel lists the number of optically thick lines (n) included in the 
least square fits, and the correlation coefficient of these fits (p). 



extinction-corrected continuum fluxes derived in the pho- 
tometric analysis in Sect. 3.2 to transform the EWs into 
absolute emission line fluxes F\. By plotting for each line: 
log(F A A 3 /(g/)) +q X u versus log(g/A) -p\h with xi and Xu 
the excitation energies of the upper and lower terms of the 
transitions, we obtained a SAC curve for each of the three 
spectra. From the 'bending point' we can obtain a lower 
limit of the Fen column density and an upper limit of the 
radius of the Fell emission region in each case. In order to 
calculate log(Ao/<?o) we need the value of vo in a thin disc, 
for which we take the thermal velocity, which is about 2.0 
x 10 cm s _1 . In the Appendix it is explained how the val- 
ues of p and q can be obtained from a comparison between 
the multiplets with common lower terms and of those with 
common upper terms respectively. 

In general the number of lines in a multiplet observed 
in our spectra was about three, but for several multiplets 
such as 27, 28, 37, 38, 49, 73 and 74 we observed more lines, 
up to 8 or 9. The resulting values of p and q show a certain 
spread related with the choice of the multiplet pairs. This 
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Figure 10. Z CMa disc radii as a function of excitation parameter 
q for the Feb. 2000 (top) and Dec. 1991 (bottom) Fen spectra. 
The insert shows the direction of varying excitation parameters, 
as discussed in the text. 



is partly caused by uncertainties in the EW measurements 
but perhaps also due to deviations from a thermal excitation 
mechanism. For the 1987 spectrum we find p = 0.9-1.13, q 
= 0.92-1.05, for the 2000 spectrum p = 0.8-1.10 and for 
the spectrum of 1991 p = 0.7-1.15, q = 0.83-1.10. Since the 
uncertainties on these numbers are ~20%, it seems justified 
to assume that there is no significant difference between the 
(p, q) values of the three spectra and that the values of p 
and q are close to 1. If the excitations of the upper and lower 
terms follow Boltzmann distributions with (excitation) tem- 
peratures defined by 5040/p for the lower terms and 5040/g 
for the upper terms of the transitions, we expect p to be 
somewhat larger than q: perhaps as large as 20% , so that 
for p = 1 (Ti =5040 K), q could be 0.8 {T u = 6300 K). 
With these choices for p and q we have calculated the SAC 
curves for the three spectra (e.g. Fig. 9 for p = q = 1) We 
note that the slope of the optically thick part of the curve 
(at the right side of the bending point), determined by least 
squares fitting, is around 0.6. Theoretical predictions of the 
SAC curve for various models of line-formation (Fricdjung 
& Muratorio 1987, Figs. 4c and 4f) indicate that such slope 
values are consistent with line-formation in winds, confined 
to a disc. 

We derived the upper limits for the disc radius and the 
lower limits of the column densities of Fe (or iV(H)) for 
a range of values of p and q, close to 1.0. Fig. 10 gives a 
survey of the resulting relative disc radii (Rd/R*) for the 
1991 spectrum, and for the 2000 spectrum for R CB = 4.2 (for 



R cs =3.1 the values of Rd/R* are smaller by a factor of 
2). Note that, although the exact size of the Fen emitting 
region in 1991 and 2000 may be uncertain by a factor 2 or 3 
(depending on the assumptions made about the excitation 
parameters), the ratio of the size of the Fell emitting region 
in 1991 and 2000 has a much smaller uncertainty, perhaps < 
50%. For (p, q) = (1, 1) when both excitation temperatures 
are 5040 K, the upper limits of Rd/R* for 1991, 2000 and 
1987 are 22, 33 and 34. This corresponds to the direction C in 
the insert of Fig. 10. For the 8542A Can (2) emission region 
we have found (previous section) that the Ca II (2) emission 
region seems to expand from Rca.n/R* = 12 to Rca.n/R* 
— 29 between 1991 and 2000. However, if we allow (p, q) to 
become lower between 1991 and 2000 the radius will increase 
less because of the rising upper term excitation temperature 
(directions B for (0.9, 0.9) and A for (0.8, 0.8)). It is also 
possible that the excitation temperature decreases between 
1991 and 2000 (direction D for (1.1 to 1.1)). In case D the 
excitation temperature has decreased by 10%, but the radius 
of the Fell disc increases by a factor of 2.0 between 1991 and 
2000. For each value of (p, q) the SAC analysis also provides 
a lower limit to the column density of Fell and (with the 
partition function and the jV(Fe)/iV(H) abundance ratio) 
a corresponding column density of Hydrogen. We used the 
abundance ratio for the solar environment. For (1, 1) wc 
find AT(H) > 7 x 10 23 cm" 2 from the spectrum of 1991. 
In the points A, B, C and D of the spectrum of Jan. /Feb. 
2000 (Fig. 10) the ratios of iV(H) with respect to the value 
of iV(H) in (1, 1) from the 1991 spectrum are 0.28, 0.52, 
1.03 and 2.06, respectively. The corresponding ratios of the 
disc radii upper limits are 0.92, 1.18, 1.52 and 1.96. In this 
example it seems that the column densities iV(H) increase 
faster than the disc radii. Since the bending point Xq of the 
SAC is independent of continuum flux, the column density 
iV(H) derived from it will not change for a different choice 
of spectral type. 

4-6.2 The forbidden Fell lines 

In the spectra of Nov/Dec. 1991 and Jan. /Feb. 2000 a small 
number of [Fell] emission lines with log(gf) > —7.9 have 
been found. The line list of Feb. 1987 (Hessman et al. 1991) 
does not contain such lines, but in the original spectra used 
by Hessman et al. such lines can also be found. The ob- 
served emission fluxes for these lines were treated in the 
same way as the permitted Fell lines and marked by filled 
plot symbols in Fig. 9. The positions of these points are 
to the left and above the bending points of the curves 
for the permitted lines. Since these [Fell] lines are opti- 
cally thin they are not formed in a dense region such as 
a disc, but probably in a extended region around the disc. 
In order to make some estimate of the extension of this 
region we applied the modified form of the SAC ordinate: 
log(V/d 2 ) = Yf - log(n /ffo) + 16.977 in which Yf is the 
normalized emission flux of the line and no is the ground 
level volume density, which can be estimated from the sur- 
face level density and radius, determined from the permitted 
lines (Baratta et al. 1998) . 

The results are presented in Table 6, which shows that 
(for p — 1.0, q — 1.0) the formation region of the [Fen] 
lines is somewhat more extended than the disc, also in the 
plane of the disc. In the outburst state of Jan. 2000 it has 
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Table 6. Radii of [Fen] emitting regions in the 1991 and 2000 
spectra for (j>,q)=(l. 0,1.0). Rd = the disc radius. Typical errors 
in these numbers arc ~ 20%. 





4385 A 


4815 A 


5158 A 


5261 A 


5333 A 


7155 A 


1991 R/R d 




2.9 


2.7 


1.8 




1.5 


2000 R/R d 


2.3 


2.0 


2.3 


1.8 


2.0 


1.2 



expanded together with the disc. For the strongest [Fell] 
lines the excitation region extends out to nearly three times 
the upper limit of the disc radius. This result appears to be 
valid for other values of p and q as well. 

4.7 The forbidden emission lines of O I, Can, Nl, 
N II and S II 

There exist at least six recent profiles of the [O i] 6300.3 A 
emission in Z CMa: the profile of March 31, 2002, a profile, 
obtained on Jan 23, 1998 (Garcia et al. 1999), one profile 
from the WHT spectrum of Dec. 1996, 4 profiles from Jan. 
14 to 17, 1997 (Chochol et al. 1998), two profiles from Dec. 
20 and 21, 1991 and one profile from Feb. 12, 1987 (Hessman 
et al. 1991). Fig. 11 shows the profiles for five dates. 

Garcia et al. (1999) were able to separate the profiles 
from the primary (Herbig BOe) and secondary (FUor) star. 
The deconvolved image in their Fig. 3 indicates that the 
contribution of the primary to the flux in the central peak 
of the [O i] emission is at least 2.5 times the contribution 
of the secondary (FUor) component of Z CMa. This is con- 
firmed by the fluxes of the recovered spectra in their Fig. 2. 
In addition they showed that the shoulder at the blue side of 
the line is due to a small (mini-)jet, with a length of ~1" in 
approximately the direction of the large-scale outflow (Po- 
etzel et al. 1989) and originating from the primary. Garcia 
ct al. suggest that this mini-jet was ejected at the time of 
the outburst of the primary in Feb. 1987. 

In our new [O i] spectra, we note a high- velocity com- 
ponent, extending blueward of the main [O i] emission peak. 
In the spectrum of 1987 we only observe a single O I peak 
without shoulder. In all subsequent [Ol] 6300 A profiles a 
shoulder is present, starting at a distance of ~ —600 km s~ 
from the peak at 6300.3 A. 

The most recent profile of the [O i] line, taken on March 
31, 2002, shows a remarkable increase in the flux of peak and 
shoulder. The appearance of these broad emission features 
suggests the sudden development of the high velocity out- 
flow component. We also note the appearance of a second 
emission peak, centered on ~ 6292 A, superimposed on the 
[O i] 6300 A line profile of Z CMa. This is similar to the sec- 
ond emission peak around 6294 A in the profile of the [O i] 
6300 A line in the spectrum of PV Cep (Fig. 1 of Corcoran 
& Ray 1997). This second peak indicates a discrete high ve- 
locity component in the outflow, perhaps due to a change in 
collimation of a jet. An alternative mechanism could be the 
emission of a discrete line or group of lines. Plausible lines 
present in the 6290-6305 A range are the lines of N II, with 
groupings of lines at 6286-6294 A and from 6298-6304 A. 

Corcoran & Ray (1997) have noted that a high velocity 
component in the [O i] 6300 A profile is a rare phenomenon 
for Herbig stars. It has only been observed for four stars: 
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Figure 11. Variations in [O i] 6300.2A line profile of Z CMa 
between Dec. 1991 and Feb. 2002. Dashed horizontal and vertical 
lines again indicate the continuum level and the rest-wavelength. 
The normalization of the profile of 2000 is somewhat uncertain 
because it is located near the edge of the echelle order. The narrow 
absorption lines seen in all spectra are due to telluric water vapour 
absorption, for which no correction was applied. 



V645 Cyg, Z CMa, PV Cep and LkHa 233. These stars are 
young, embedded stars which are associated with stellar jets 
and molecular outflows and probably accreting rapidly. The 
jets form close to the stellar surface and are ejected along 
the rotation axis of the star. The collimation of such jets 
within 150-200 AU has been successfully modeled with a 
stellar or disc-generated wind in a rapidly rotating magne- 
tosphere around such stars. The fact that such high velocity 
jets are not observed for the majority of the Herbig stars 
may be caused by differences in the evolution of the colli- 
mating magnetic field. 

We have also inspected the [Call] 7291.46 A and 7329 A 
profiles in the Tautenburg spectrum and in the spectra of 
Nov. /Dec. 1991 and Dec. 1996. It is remarkable that they are 
similar to the central peak of the [O i] 6300 A line, but they 
clearly lack the extended blue shoulder. The same seems to 
hold for the [O i] 6363 A, the [S II] 6731 A, the [N II] 6583 A 
and the [Ol] 5577 A profiles, but this is more difficult to 
verify since the flux in the peak of these lines is three times 
lower than that of the [O i] 6300 A line. It is not completely 
clear whether the central peaks of the [N n] and [S n] lines 
are mainly due to the primary BOe star. However, it seems 
convincing that the total flux is dominated by the primary 
since the [N n] and [S n] fluxes seem confined to the jet direc- 
tion, which runs through the primary (Poetzel et al. 1989; 
Garcia et al. 1999). The weakness of the [N n] and [S n] lines 
does not allow us to obtain reliable profiles. If most of the 
forbidden lines originate in the neighbourhood of the pri- 
mary, we can try to use the flux ratios to obtain a rough 
estimate of the conditions in the far outer atmosphere of 
the BOe star. Hamann (1994) calculated ratios of various 
forbidden line fluxes with respect to the line flux of [O i] 
6300 A under the assumptions of coronal ionization equilib- 
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Table 7. n e (cm' 3 ) for T e 
numbers are ~ 30%. 



=13,000 K. Typical errors in these 
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Figure 12. Fragments of the Z CMa IRTF spectrum of Dec. 
2001, showing the lines of Bra (4.05 /tm), Pf7 (3.74 /im), Br7 
(2.67 /im) and the CO band-head around 2.3-2.4 fim. 



rium and solar abundances. This allows us to estimate T e 
and n e from the flux ratios of [Nil] (6583 A), [Sir] (6731 A), 
[Can] (7291 A) and [Oi] (5577 A) to [Oi] (6300 A). We have 
used the observed EWs and intrinsic continuum fluxes after 
circumstellar extinction corrections with R cs = 4.2. From 
the diagrams of Hamann we then obtain for T e — 13,000 K 
the values of n e in the formation region of the various lines 
(see Table 7). 

For T e = 15,000 K we find the same trends: (a) for 
a fixed T e the [S n] line comes from the lowest density 
and therefore at the largest distance from the stellar en- 
velope, (b) the densities decrease when the Be star becomes 
brighter. This suggests an expansion of the outer region 
where the forbidden lines are formed. 



4.8 The IRTF spectrum 

In the IRTF spectrum of Dec. 2001 we observe the emis- 
sion lines Bra (4.05 /mi), Pfy (3.74 fim), Br7 (2.16 /im). 



Table 8. Lines in the IRTF spectrum 



A [/im] 


Ident. 


EW [A] 


FWHM [km s" 1 ] 


Dout [km s : ] 


4.05 


Bra 


-7.75 


370 


-1600 


3.74 


Pfy 


-6.37 


545 


-1670 


2.37 


CO (2-0) 


4.28 


328 




2.32 


CO (2-0) 


2.31 


328 




2.17 


Br7 


-5.00 


375 


-1500 



He I (2.06 pm), Br<5 (1.94 /mi) and absorption in the CO 
v — 2-0 overtone bands around 2.3 /im. Br/3 is in the ob- 
served wavelength-range, but it is obscured by terrestrial 
atmospheric bands. Fig. 12 shows the profiles of Bra, Pfy, 
Br7 and the CO band. The EWs are given in Table 8. The 
CO absorption most likely arises in the FUor component of 
Z CMa (Hartmann et al. 1989). The infrared Hi emission 
lines appear to have similar outflow velocities as their opti- 
cal counterparts (Fig. 1); they may be related to the inner 
accretion disk around the Herbig B0 primary. Muzerolle et 
al. (1998b) have derived empirical relationships between the 
luminosity of Br7 and the accretion rate of T Tauri stars. 
Our continuum flux at 2.2 /mi corresponds to K — 4.03 m . 
We therefore have a Br7 emission line flux of 6.0 x 10~ 12 
erg cm" 2 s" 1 or logi ? (Br7) = —11.2. From Fig. 3 in Muze- 
rolle's paper we then find by extrapolation that the accretion 
rate of Z CMa in Dec. 2001 should be ~2 x 10" 5 M Q yr" 1 . 

4.9 Summary 

In the preceding sections we have compared the spectra of 
the 'outburst' phases in Feb. 1987 and Jan. /Feb. 2000 with 
those of the 'quiescent' phases Nov. /Dec. 1991 and Dec. 
1996. The most remarkable distinction between these phases 
is the appearance of deep absorption lines of He I at 7065, 
6678 and 5876 A and of O I at 7773 A and 6156 A in the out- 
burst spectra. We interpret this appearance as a clear indica- 
tion that the 'outburst' occurs in the primary Be component 
of Z CMa. 

The strong similarity between the spectra of Z CMa 
with those of the class of Herbig Be stars suggests a similar 
interpretation: the emission lines are formed in a circum- 
stellar disc. Because of the direction of the huge bi-polar jet 
(large radial velocities in the knots), the association of this 
outflow with the primary star (Garcia et al. 1999) and the 
observed polarization of the lines (Whitney et al. 1993), we 
assume that we are seeing the disc under an angle of at most 
40 degrees. 

Following Hamann & Persson (1992b) we have used the 
fluxes of the strong Can (2) triplet to estimate the extent 
of the corresponding excitation region. We found that the 
outer radius of this region expands from 10 R* in 1991 and 
1996 to ~ 29 R* during the outburst in Jan. /Feb. 2000. The 
line profiles of Call (2) may indicate that during the faintest 
phase (Dec. 1996) a small contribution of Call emission from 
the FUor component may be present. The analysis of the 
Fell emission multiplets shows a similar expansion of the 
Fen region. This expansion is probably somewhat smaller 
than that of Call: 23 R* in 1991 to 32 R* in 2000, but 37 
R* in Feb. 1987. [Fen] emission extends up to twice this 
distance. Other forbidden lines in the spectra such as [Sll], 
[N n] and [Ca n] correspond to very low values of n e and are 
formed much further away, probably in the bipolar outflow 
(Poetzel et al. 1989). 

The similarity of the Balmer and Paschen outflow ve- 
locity pattern with those of the Call (2) and Fell (42) lines 
(both groups have several outflow components with veloc- 
ities up to 450-600 km s" 1 ) suggests that these lines are 
formed in the same region. The emission components of Ca II 
(8542 A), Pll (8863 A) and Oi (8446 A) are very similar in 
profile. Especially in the 1991 spectrum, where the blue ab- 
sorption wings are less conspicuous, the profiles are very 
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similar to those of the T Tauri stars BP Tau and UY Aur. 
Muzerolle et al. (1998) have modeled these BP Tau pro- 
files with a nonrotating, axisymmetric dipole field for mag- 
netic accretion and gave rough relationships between the 
extinction-free line fluxes of the Call (8542 A), Ol (8446 A) 
and Pll on one hand and the accretion rate derived from 
the UV/blue excess continuum (Gullbring et al. 1998) on 
the other. We have used these relationships to derive for 
Z CMa accretion rates of logM (M yr -1 ) = -4.4 for the 
Ol and Pll and —5.3 for the Call emission lines by extrap- 
olation of Fig. 9 of Muzerolle et al. (1998). From the flux 
of Br7 and Fig. 3 in the paper of Muzerolle et al. (1998b) 
we obtain by extrapolation that logM (Mq yr _1 ) = —4.7. 
It therefore seems possible that in Z CMa these lines show 
inflow (magnetospheric accretion) as well as outflow. 

Because of the changes in line emission flux and outflow 
velocities we expect that both flows are different in the var- 
ious phases of Z CMa. The He I absorption lines have high 
excitation energies and must therefore be formed close to 
the star while the Ol absorption lines probably have their 
origin in an extended envelope above and below the disc 
plane. Most of the emission lines of O I, Can, Fen, Fel, 
Cm, Mgn and the Paschen lines form in the disc. NalD, 
K I and Ca II K are formed further out in the envelope. The 
CallK absorption pattern shows similarities with those of 
the NalD lines. These lines may receive contributions from 
the envelope of the FUor component as well. However, it 
seems highly probable that the variations in these profiles 
also reflect changes in the outer envelope of the primary Be 
star during the outburst phases. 

Two new phenomena have been observed in connection 
with the outburst of Jan. 2000: the appearance of a strong 
Fen 9997 A emission line and of an extended blue shoulder 
in the [O i] 6300 A emission profile in 2002. The origin of the 
first emission is not yet clear. It is possible that Lya fluo- 
rescence is the mechanism for its production. The emission 
profiles of the [Ol] 6300 A and the [Can] 7291, 7323 A dou- 
blet profiles show a gradual growth in EW since Dec. 1987. 
The blue shoulder in the [O i] 6300 A profile, detected two 
years after the outburst of Jan. 2000 is possibly caused by 
a sudden change in the collimation process of the mini-jet, 
which could have been ejected in Feb. 1987. 



5 DISCUSSION 

The spectrum of Z CMa in its normal state has been satis- 
factorily explained by Hartmann et al. (1989), and Welty et 
al. (1992) in terms of that of a rotating, optically thick, ac- 
cretion disc. From a total luminosity of 3000 Lq Hartmann 
ct al. (1989) derived for «sini ~ 100 km s _1 an inner disc 
radius of 16 R and an accretion rate of ~ 10~ 3 Mq yr _1 
to a central star of 1-3 Mq and radius ~ 9-16 R©. How- 
ever, during the outburst of Feb. 1987 (Hessman et al. 1991) 
the spectrum was completely different and the observers of 
this high state could not explain it in terms of the same 
or an adapted model. It were the spectropolarimetric ob- 
servations of Whitney et al. (1993) that showed that the 
usually small contribution of the infrared primary star to 
the visual continuum could rise to ~66% of the total contin- 
uum flux during the outburst of 1987. Whitney et al. (1993) 
have noted that in this high state the emission spectrum of 



Z CMa is very similar to that of the Herbig Be star MWC 
1080. The agreement is not limited to the emission lines but 
also to the profile of the NalD line (Finkenzeller & Mundt 
1984). Similar to Z CMa, the He I lines at 6687 and 5876 A 
of MWC 1080 are in absorption (Yoshida et al. 1992). 

We also note a strong spectral similarity with the lu- 
minous young object V645 Cyg (Hamann & Persson 1989), 
which has been classified as 08. Especially the broad Fen 
9997 A emission line and the broad Ol (7774 A) absorption 
line together with a relatively weak O I (8446 A) emission is 
very similar to the high state of Z CMa. This indicates an 
important role of Lya and Ly/3 emission in both stars. A dif- 
ference is that the He I (6678 A) P-Cygni profile of V645 Cyg 
has a more prominent emission component than its counter- 
part in Z CMa. This may point to a somewhat hotter inner 
disc. Another feature which both objects have in common 
is the high velocity component of the [O i] 6300 A emission 
profile (see below). 

The spectrum of Z CMa has also been compared with 
the rich emission line spectrum of V380 Ori. However, the 
spectrum of V380 Ori (Rossi et al. 1999) shows some marked 
differences with those of Z CMa and MWC 1080. Its He I 
lines at 6678 and 5876 A are generally in emission, but occa- 
sionally in absorption. It has been suggested by Rossi et al. 
and others that this may indicate the presence of a chromo- 
sphere or a hot accretion spot close to the stellar surface. Its 
O i (7773 A) line is also in emission, which is not unusual for 
various early A-type Herbig stars, such as AB Aur (AOe), 
HD 163296 (A2e) and BD +61° 154 (B8e) (Felenbok et al. 
1988). An independent argument for a classification of V380 
Ori as a B9-A0e star is the cut-off of the ultraviolet contin- 
uum near 1300 A in the low resolution IUE spectra of V380 
Ori (Rossi et al. 1999). 

The spectral type of the Z CMa primary component 
was found to be close to BOe (Sect. 4.2). From V = 6.05 m 
and a distance of 1050 pc we find an absolute magnitude 
Mv = — 4.1 m , which may indicate a classification Bl-2 III. 
In Sect. 3.2 we have derived the luminosities of the primary 
component for Nov. 21, 1991 (the date of the polarization 
spectra of Whitney et al. 1993) and Feb. 20, 1987 (the 'high 
state' of Z CMa) . The results for three assumptions for R cs 
(the ratio of total to specific absorption of the circumstellar 
envelope): 3.1, 4.2 and 6.0 are given in Table 2. 

Is may be interesting to compare the temperatures and 
the de-reddened luminosities of this table with recent calcu- 
lations of the pre-main sequence evolution of young massive 
stars by Behrend & Maeder (2001). The comparison shows 
that in the high state for R CB = 3.1 and R CB = 4.2 the emis- 
sion line star can be identified as a pre-main sequence star, 
since then it can be modelled as a 16 Mq B0 III star (T e ff 
~ 31,600 K) on the birthline with an age of 3 x 10 5 yr and 
a luminosity of 47 L and radius 7.6 R . For R CB = 6.0 the 
luminosity is too high and the position of the star is above 
the birthline and, in order to be observable, it should be a 
post main sequence star. From the post main sequence evo- 
lutionary calculations of Schaller et al. (1992) we find that 
it could be 20 M Q B0 III star with an age of ~ 6.6 Myr. 
This is close to the estimated ages of the nearby BOe stars 
HD 53367 and HD 53755 (paper II). On the other hand the 
parameters of the inner star of the FUor component, esti- 
mated by Hartmann et al. (1989), correspond to a 3 M star 
with a radius of 7 Rq, an effective temperature of 5000 K, 
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a luminosity of 28 L0 and an age of 3 x 10 5 yr. Its position 
in the HR diagram is just below the birthline. In this case 
both components of Z CMa have the same age and could 
have been born together. 

Recently Yorke & Sonnhaltcr (2002) have published 
evolutionary models for star formation in which radiative 
acceleration of dust grains has been taken into account in 
the radiative transfer. The results of the evolution depend 
on the opacity law of the dust. For a frequency-dependent 
opacity law due to a mixture of two dust components with 
specified size distributions the evolution leads to the for- 
mation of massive stars via disc accretion with a massive 
bipolar jet. For a 'gray' opacity law, however, the evolution 
proceeds without bipolar jet and results in a thin, disc-like 
object. Since in these models the birthline is raised by a 
factor ten (for BO stars) with respect to the previous mod- 
els, they can easily interpret the Z CMa primary luminosity 
of Feb. 1987 with R cs up to 6 as the luminosity of a pre- 
main sequence star. The paper shows 2-dimensional models 
of the immediate neighbourhood of the stars for a wide grid 
of masses, ages and luminosities. 

For our case of 3.1 x 10 5 L (Table 2, R cs = 6) we are 
not far from a model of a 38 Mq star, close to the birthline 
with an age of 4 x 10 4 Myr (model F120, 11c). For this age 
the model predicts a mass accretion rate smaller than 10~ 4 
Mq yr -1 , which is consistent with our estimate from the Br7 
flux in Sect. 4.8. Similar to Z CMa, the 2-dimensional model 
shows an indication of an accretion disc and of a strong bi- 
polar outflow. The dynamical time scale of 2-3 x 10 4 years 
for the 'slow' knots (Poetzel et al. 1989) in the bipolar jet 
fits within the model age of this BOe star. 

For ages comparable with those of the BO star models 
the FUor component could correspond to a 'gray' opacity 
model such as G60 of Yorke & Sonnhalter (2002) but for a 
lower luminosity and mass than calculated in their Table 4. 
Such a model does not develop a bipolar outflow, but evolves 
to a disc-like structure, perhaps a FUor. In this case we may 
have a simultaneous birth of the Z CMa components. This 
is in agreement with the polarization results of Fischer et al. 
(1998), who came to the conclusion that both components 
must have had a common formation history. However, the 
models by Yorke & Sonnhalter imply a difference in grain 
properties for the regions in which the B0 and the FUor com- 
ponent are formed. Because the present rotation velocity of 
the B0 component and its previous evolution are not known 
we cannot be certain of the evolutionary age, but so far it 
seems that the age estimates of Z CMa with both models 
are much lower than those of the other massive stars in the 
CMa Rl complex, which are closer to the main sequence 
(paper I, II). How can the low age of Z CMa be fitted in 
the evolutionary history of CMa Rl? The formation must 
be related to fragmentation in the high density in the dark 
cloud S296 where Z CMa is situated, but so far it remains 
unclear which process initiated the contraction in the cloud 
which started the recent star formation. 

Apart from the long-term evolution of Z CMa there re- 
main a number of questions concerning the recent behaviour 
of the two components: 

(1) The cause of the photometric variability: At least 
two types of variabilities have been observed: (a) irregu- 
lar low amplitude variations, which seem on the average 
to follow density variations in the circumstellar dust. V.S. 



Shevchenko made an extensive search of periodicities in the 
UBVR photometry, collected by the ROTOR program dur- 
ing the years 1981-1998. A period of 393 days was found 
with an amplitude of 0.3 m in U - B and 0.08 m in B - V. 
The underlying cause of such a periodicity is not yet known 
but pulsational instability may be a possibility, (b) Strong 
and fast rises in brightness (outbursts) around 1968 (AV ~ 
0.7 m ), 1987 (AV ~ 0.6 m ) and 2000 (AV - 1.2 m ). We have 
now some indications that the second outburst has caused 
the ejection of the mini-jet emitting [0 1] 6300 A (Garcia et 
al. 1999), similar to those observed in other young emission- 
line objects such as PV Cep, V645 Cyg and LkHa 233 (Cor- 
coran & Ray 1997). The third outburst (in Jan. 2000) prob- 
ably caused a change in the [O 1] 6300 A profile (see Fig. 11) 
which could be due to a change in the collimated outflow. 
The outbursts therefore could be related to changes in the 
upper layers of the stellar interior, which entail changes in 
the magnetic field configuration. Such changes could have 
an internal origin or be induced by an origin from outside 
the star, e.g. by tidal interaction during approaching compo- 
nent stars. The latter possibility (Bonnell & Bastien 1992) 
can be ruled out by the following argument: The projected 
distance between the components of Z CMa is ~ 0.1" which 
at a distance of 1050 pc corresponds to ~105 AU. This is 
seen in projection and therefore the minimum distance of 
the components. A minimum for the period can then be es- 
timated from Keplers law. For Mi = 38 M Q and M 2 = 3 M , 
we find a minimum period of 168 yrs. This is too large for an 
observed recurrence of outbursts at a timescale of once per 
10-12 yrs. The remaining possibilities are that the primary 
star has a, as of yet undetected, close binary component 
or that the outbursts are induced by variations in the star 
itself. 

(2) According to Garcia et al. (1999) the direction of 
the mini-jet coincides with that of the large bi-polar outflow 
of Z CMa and the outflow is driven by the massive B0 com- 
ponent. Also the evolutionary models of Yorke & Sonnhalter 
(2002) predict large bi-polar outflows during the formation 
of massive stars. However, Velazquez & Rodriguez (2001) 
claim that the large bipolar outflow coincides with the radio- 
jet through the FUor component. If this is true, the question 
remains how the cavity in the cocoon is formed, through 
which we observe the optical spectrum of the B-star. One 
possibility is that it has been formed in the past when the 
infrared source approached the FUor. Somehow the direc- 
tion of the bipolar outflow is then conserved in the infrared 
source and transferred to the mini-jet. 
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APPENDIX A: THE SELF-ABSORPTION 
CURVE (SAC) METHOD 

The observation of many multiplets of emission lines of Fe II, 
Ti II and Cr II (Table 4) allows us to derive some information 
concerning the formation regions of these lines with the help 
of an emission curve of growth method. The method is based 
on the analysis of the self-absorption curve (SAC) which 
was developed by Friedjung & Muratorio (1987) and has 
been successfully applied to various Fe II emission line spec- 
tra (e.g. Muratorio et al. (1992) for the case of KQ Pup)§. 
A manual for the use of the method has been written by 
Baratta et al. (1998). 

The SAC curve describes the relation between the line 



§ An independent derivation of the SAC method has been given 
by Kastncr (1999) 
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strength and optical depth of emission lines in a multi- 
plet. This relation can be empirically determined by plotting 
log (_F\ A 3 /gf) (where F\ is the normalised absolute line flux 
and gf is the oscillator strength of the line) versus \og(gf) 
(which is the optical depth at line-center up to an addi- 
tional constant) for the various lines of each multiplet. In 
these empirical plots the lines through the points of each 
multiplet should be parallel to each other and to the still 
unknown theoretical SAC function Q(t), which should be 
zero for log r < and should should decrease approximately 
linearly for logr > 1 . The emission lines in the horizontal 
part of the SAC (with logr < 1) are optically thin, whereas 
those in the decreasing part are increasingly optically thick. 
For the various Fen emission lines the oscillator strength 
and excitation energies have been taken from the tables of 
Kurucz (1981). 

The various curves for multiplets with common upper 
terms ^ will be parallel and shifted by pxi in horizontal direc- 
tion with respect to each other, whereas those with common 
lower terms 1 1 will be shifted by q\ u in vertical direction with 
respect to each other. Here \i an d Xu are the excitation po- 
tentials of the lower and upper terms in the transitions and 
p and q are the values of 5040/T; and 5040/T„, where it is 
assumed that the levels have Boltzmann-type populations 
with excitation temperatures Ti and T u for the lower and 
upper terms respectively. Since we know the excitation po- 
tentials, we can determine the values of p and q from the 
horizontal and vertical shifts and therefore the Boltzmann 
temperatures. 

By plotting Y = \og{F\\ i /gf) + q\u versus X = 
l°g(ff/^) — PXi f° r the lines of each multiplet, we can obtain 
one common empirical SAC curve for the Fen lines. The 
value Xo on the X-axis, at which the empirical SAC curve 
starts to decrease for increasing X, gives us log(Nofgo) up 
to constant log« + 1.576 (in cgs units), where v is the line 
broadening velocity and No is the ground state column den- 
sity of Fe II from which the total Fe and H column densities 
can be estimated. From the normalized emission line flux Yo 
in this point Xo, where the empirical SAC curve starts to 
decrease, it is then possible to estimate the extension of the 
surface area S of the emission region (perpendicular to the 
line of sight) using log(S/d 2 ) =Y -log(N Q / g )-w + 16. 977), 
where d is the distance to the source (here 1050 pc) and 
w — —(p — q) x (xi + X«)/2- An example of the empirical 
SAC curves of Fe II emission lines of three spectra of Z CMa 
(observed in Feb. 1987, Jan. 2000 and Nov. 1991) is given 
in Fig. 9. 

For spectra with only optically thick lines, one can only 
obtain the decreasing part of the SAC curve, so that the 
values of Xo and Yo have to be estimated. The method then 
only gives a lower limit to the column density and an upper 
limit to the extension of the emission line region. On the 



other hand the method can be applied to optical thin for- 
bidden lines, to derive a characteristic volume (or radius) of 
their emission region. Finally the slope of the decreasing part 
of the theoretical SAC curve has been shown to depend on 
the dynamical conditions in the emission region (Friedjung 
& Muratorio 1987). In our case of Z CMa the slope is close 
to —0.6, which may indicate that the Fell lines are formed 
in a disc wind (Fig. 4f in Friedjung & Muratorio 1987). 

This paper has been produced using the Royal Astronomical 
Society /Blackwell Science IATpX style file. 



^ In the visual part of the spectrum Fe II multiplets with com- 
mon upper terms are: 27, 38, 48 and 73 ( 4 P( 5 D) 4 D), 28, 32, 37, 
49 and 55 ( 4 P( 5 P) 4 P), 26,42 and 36 ( 4 P( 5 P) 6 P), 29 and 74 
( 4 P( 5 D) 4 P), and 41 and 46 ( 4 P( 5 D) 6 P). 

II In the visual part of the spectrum Fell multiplets with com- 
mon lower terms are: 27, 28 and 29 (D 6 S ( 3 P) 4 P) , 36, 37 and 38 
(D 6 SC i F) i F), 40, 41, 42 and 43 (P 5 5 2 ( 7 5) 6 5), 46, 48 and 49 
(D 6 5( 3 G) 4 G) and 72, 73 and 74 (D 6 S ( 3 D) 4 D) . 
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